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Preface 


The purpose of this thesis has been to present a self- 


contained text explaining the fundamentals of the high altitude 


electromagnetic pulse in a manner which would provide a basic 


understanding of the phenomenon for the first year graduate 


student of Engineering. 
I would like to emphasize the fact 


that it has not been the purpose of this thesis to provide 


"state of the art" threat levels for the fields of the 


electromagnetic pulse. 
Rather, the calculations concerning 


these levels which are included in this report are meant 


only to show order of magnitude agreement of the model 


developed here with accepted theories. 
Although some of 


the results of derivations in the report involve advanced 


electromagnetic theory, 
I 
have attempted to develop these 


derivations in a manner that requires only a familiarity 


with integral and vector calculus, Maxwell's equations, and 


a few basic precepts of modern physics. 


I would like to express my sincere appreciation to 


Major Carl T. Case for the many hours of guidance and advice 


which he devoted to this effort. 
I would also like to thank 


Dr. Charles J. Bridgman for his encouragement and interest 


in the project. 
Finally, 
I would like to acknowledge my 


wife's patience and support, without which I could not have 


completed this thesis. 
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Abstract 


This thesis constitutes a microscopic approach to the 


electromagnetic pulse 
(EMP) 
caused by a high altitude nuclear 


weapon burst. 
This approach involves the detailed calculations 


for the time-dependent electromagnetic fields of a single 


relativistic, 
Compton electron created in the upper atmosphere 


by the gamma radiation from the weapon detonation. 
The effect 


of all such electrons created is superimposed to obtain the 


net effect or the EMP measured by an observer. 
The basi<=: 


theory developed ignores the conductivity which builds up 


in the surrounding air. 
This is shown to be a relatively 


good assumption for weapons with gamma radiation yields less 


than a kiloton. 
A brief section is included to indicate the 


derivations necessary for the EMP in a conducting media. 
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THEORY OF THE HIGH ALTITUDE 


ELECTROMAGNETIC PULSE 


I. 
Introduction 


With the advent of nuclear weapons several previously 


unsuspected phenomena emerged as possible system damage 


mechanisms. 
One such phenomenon is the interaction of the 


radiation from a nuclear weapon with the environment in 


such a mapner as to cause the motion of electric charges. 


The motion of these charges in turn causes a pulse of elec- 


tric and magnetic fields which can be very intense and can 


be radiated for great distances. 
The nature of these fields 


varies according to the altitude at which the weapon is 


detonated. 
This thesis will concern only those fields 


associated with weapon detonations which occur at exoatmos- 


pheric altitudes. 
With the development of weapon systems 


involving highly sensitive electronics, the electric and 


magnetic fields of the electromagnetic pulse 
(EMF) 
asso- 


ciated with high altitude bursts constitute a phenomenon 


which should be fully understood by defense planners. 


Definition 


The EMP can be defined as a long range nuclear weapon 


effect wherein a high intensity, radio-frequency electro- 


magnetic pulse is generated 
(Ref 9, chap X). 
There are 


several other similar nuclear effects which should not be 


confused with EMP. 
These include the Argus effect, the 


blackout phenomenon, and transient effect on electronics 
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i 
(TREE). 
The Argus effect is the phenomenon whereby the 


electrons produced by a nuclear burst are trapped in the 


earth~s magnetic field thereby creating a potential radia- 


tion hazard for personnel and electronic equipment 
(Ref 2, 


pg 1). 
The blackout phenomenon relates to the ionization 


produced by a nuclear detonation in or above the atmosphere 


and the subsequent radio frequency interference 
(BLACKOUT) 


resulting from the attenuating, reflecting, cluttering and 


scatteri~g of radar and radio signals 
(Ref 2, pg 1). 
TFEE 


refers to the direct effects of radiation on electronic 


systems and ·components 
(see also Ref 14, pp 1-3). 


History 


Although it has been said that Enrico Fermi predicted 


the EMP before the Trinity explosion on 16 July 1945 
(Ref 14, 


pg 4), the high frequency component of this signal associated 


with high altitude weapon detonations did not become a seri- 


ously studied problem until the early 1960's. 
There are 


several reasons for the recent interest in this problem. 
One 


of these is the fact that the electronics in weapons systems 


have become more sensitive to extraneous electrical fields. 


Furthermore, it has been shown that the high altitude EMP 


can serve as a "signature" for a weapon detonation, thereby 


allowing for greater control of the ban on atmospheric test- 


ing. 
Also, there has been a greater understanding of the 


impact of the effect on systems previously thought hardened 


by defense planners. 


Some of the first serious studies of the problem of the 
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high altitude EMP began in 1961 with theoretical studies by 


w. J. 
Karzas and R. Latter 
( Ref 12). 
These studies pro- 


vided theoretical models and mechanisms for describing what 


was then the unconfirmed high frequency component of the high 


altitude EMP. 
The existence of this component was confirmed 


in the 1962 Pacific tests in which "a component of the 


'radio flash' 
which depends on the magnetic bearing of the 


nuclear detonation from the observation point" was observed 


(Ref 20, pg 3). 
From the early 1960 us until the present, 
• 
the EMP, and more specifically the high altitude, high 


frequency EMP, has received a great deal of theoretical 


and experimental study. 


High Altitude EMP Source 


The EMP which results from a high altitude weapons burst 


is generated by electrons which are created in the upper at- 


mosphere. 
These electrons are created by energetic photons 


(y-rays) 
from the weapon which undergo collisions with the 


molecules of the air, ejecting atomic electrons. 
These 


scattereo electrons are very energetic and radiate by means 


which will be described in a later chapter. 


It should be noted that electromagnetic signals can 


arise from a nuclear burst from sources other than that 


mentioned above. 
An example of one such signal is the case 


signal which is associated with assymetries in the weapon 


casing. 
Another involves the electromagnetic fields gen- 


erated within a system by the interaction of the weapon 


radiation with the materials of the system itself 
(System- 
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. 
. 
Generated EMP or SGEMP). 
A third example concerns the 


electrostatic fields in the near vicinity of the burst 


caused by charge separation due to ionization 
(Close-In 


EMP). 
The relative threats of these other effects are small 


as compared to the high altitude EMP. 


Microscopic Versus Macroscopic 


When calculating field strengths associated with the 


high altitude EMP, the general approach is "macroscopic" 


in nature. 
That is, one solves Maxwell's equations for 


electromagnetism using expressions for the "current" of 


electrons created in the upper atmosphere. 
Calculations 


are made based on the "current density" or the amount of 


charge passing through a unit area in a unit time interval. 


If one can find an expression for the current density, then, 


in principle, one can compute. the electric and magnetic 


fields associated with this current. 


Although the approach described above has been used 


with much success to obtain actual values for field strengths 


associated with the EMP, it is not an enlightening approach 


for one who is interested in the physical processes which are 


involved in the generation of the EMP. 
There is a consid- 


erable amount of physics "buried" in this macroscopic approach 


to the EMP. 
For this reason, it is the purpose of this thesis 


to develop a microscopic model of the high altitude EMP. 


Rather than calculating the EMP field strengths based on a 


current density, the fields for a single, accelerated electron 


have been calculated in great detail. 
The effect of all the 
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electrons created are then superimposed to obtain the over- 


all or net field strength of the generated EMP. 
It is hoped 


that-this microscopic approach will enable the serious stu- 


dent of the EMP to have a more basic and satisfying under- 


standing of the high altitude EMP. 


Overview 


The first chronological step in the generation of the 


high altitude EMP is the production of gamma rays resulting 


from a nuclear burst. 
These gamma rays·travel through 
t~e 


upper atmosphere depositing energy through Compton collisions 


thereby creating energetic electrons 
(See Figure 1). 
In 


Chapter II, which deals with the nature of the weapon gamma 


radiation and its transport through the atmosphere, it is 


shown that the number density of Compton electrons, Npri ' at 


a distance R from a burst occurring at an altitude of HOB 


is given by 


(11-17) 
N 
- (R,t) 
= 
A (R) 
prJ.. 


y 


X 
{ 
Y 
Ey 


tf f(t'- R ) dt'} 
o 
C 


where A is the Compton macroscopic 
cross section, 
~ 
is the 
o 
total macroscopic cross section at sea level, 
Z 
is the 
0 
scale height, 
y 
is the total gamma yield, 
E 
is the average 
y 
y 
gamma ray energy, f(t) is the time dependence of the yield, 


and an exponential atmosphere has been assumed. 


The Compton electrons are turned by the earth's 


magnetic field which causes them to radiate. 
It is shown 
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in Chapter III that, 
for a given observer, the radiated- 


fields which he measures result from only those electrons 


moving within a small cross-sectional area centered about 


his line-of-sight with the weapon burst 
(see Figure 1). 


Further, it is shown that, for the periods when its radiation 


is significant, the electron has essentially a constant speed. 


The power radiated by such an electron is primarily in the 


forward direction 
(see Figures 15 and 16) as is indicated 


by the following formula 
(developed in Chapter III) 
for the 
dP 
.. 
power radiated per solid angle 
(dr2C) : 


dP 
. 
2 
sin 2l/J cos 2 
(111-40) 
c 
_ 1 
(~) 
1 
{l - 
n 
} 
dr2 
- 
C£ 
41T 
(l-S cos l/J) 3 
y2(1-S cos l/J) 2 
V 


where C, 
£v' 
e are the speed of light, free space per- 


mittivity, 
and electron charge respectively, 
S lf the 


electron's relative speed, S is its relative acceleration, 


and 
l/J and n are angles defined in Figure 14. 
If A(R) 


represents the cross-sectional area which encloses electrons 


contributing to the EMP measured by an observer located at a 


distance X from the burst, then the time dependent fields 


which he observes are given by: 


a 
(III-46) E(X,t) = Yy fdR A(R)A(R)F(R) 
Ey 
0 
(X-R) 


t- (X-R) 
C 
J 
dT"f(T")E (t-T") 
e 


where F(R) is the R dependence given by the exponential term 


divided by 
41TR2 in 
(11-17) 
above and E (t) is the general 
e 
time dependence of the field strengths radiated by a single 


electron and measured by the observer. 
This equation is 


developed in Chapter III. 


Determining the size of A(R) is one· of the key cal- 
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culations in this approach. 
Chapter 'IV illustrates 
ho~ to 


obtain this variable.and also gives some results of the EMP 


which an observer would measure. 
These calculations were 


performed assuming that the conductivity of the air was 


zero. 
The range of validity for this assumption is checked 


in Chapter V and it is indicated how one might proceed 


when air conductivity is significant. 
The equations above 


and the developments which follow are relativistically 


correct. 


• 
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~~EAPON BURST 


i ELECTR01~S 
TURNED 
BY 
" ~~HTE' S 
r.AGU~TIC FIELD 


\ 


. ! 


- 
'¥t=2 


RADIATION PATT~I 
FOR S n~GLE ELECTROIJ i 


II 
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Figure 1. 
Generation of the High Altitude EMP. 
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II. 
Review of Nuclear Weapon Physics and Effects 


Since the EMP depends upon the radiation emitted from 


a nuclear weapon, it is the purpose of this chapter, to pre- 


sent a brief discussion of the physics involved in the creation 


and transport of the radiation produced from a nuclear weapon 


detonation. 


The chapter will encompass three main areas relating to 


weapon radiation. 
The first area will involve the origin 


and nature of the radiation that is proquced within the 
. 


initial weapon burst. 
Secondly, the transport of this 


radiation through the surrounding atmosphere will be con- 


sidered. 
Finally, the interaction of the radiation with 


the atmosphere to create energetic electrons 
(Compton effect) 


will be investigated. 
Key derivations are included in the 


Appendices and are referenced throughout the main body of 


the text. 


Origin and Nature of Radiation 


With regard to the high altitude EMP the important 


portion of the weapon output is the radiation in the form 


of gamma rays and x-rays. 
For weapon detonations occurring 


below 100 or so kilometers, the primary source of the EMP 


involves only the electrons created by the gamrna rays. 
One' 


reason for this is that at these altitudes the range of the 


electrons created by the x-rays is too short to cause an 


appreciable current. 
Furthermore, for detonations at 


these altitudes, the x-rays are absorbed in a relatively 


small region about the burst point. 
This region has been 
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made conducting by ionization from the earlier emission of 


gamma rays and hence the weak radiation from the electrons 


created by the x-rays is effectively damped out 
(Ref 11, 


pp 40-5). 
Since the basic physics for the creation of the 


EMP is essentially the same for x-rays and gamma rays, only 


the EMP generated by gamma rays will be considered in this 


thesis. 


There are two main characteristics of the gamma radi- 


ation output of the weapon which must be considered with 
• 
respect to the EMP. 
First it is necessary to know the 


energy of the gamma rays and second, one must know the time 


dependence of the radiation yield of the weapon. 
The energy 


of the gamma rays produced in a typical theoretical fission 


device range from under 1.0 MeV up to around 5.0 MeV, with 


1 
an average energy of approximately 1.5 MeV. 
For simplicity, 


in any calculations which follow it shall be assumed that the 


spectrum of gamma rays emitted from a weapon is monoenergetic 


such that all gamma rays have an energy of 1.5 MeV. 
The time 


dependence of the radiation yield of a hypothetical weapon is 


characterized by a sharp rise followed by a decay both of 


which are generally described as exponential in nature. 
The 


physical reasons for the time dependence of the radiation 


output is beyond the scope of this study, however, the in- 


terested reader is referred to The Physics of Nuclear 


lInformation concerning the average energy was obtained 
in private correspondence with Dr. 
C. J. Bridgman, Physics 
Dept., AFIT. 
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Explosives which is a set of class notes by 
Dr~ c. J. 


Bridgman 
(Ref 1) . 


.There have been several attempts to analytically 


describe the time dependence of weapon yield. 
One such 


attempt which is an empirical model for the time depend- 


ence presents the following formula. 


(II-I) 
f (t) 


(a+S)exp{a(t-t )} 
o 
= N(B+a exp{ (a+B) (t-t )}) 
o 
(Ref 18) 


where a and 
B are weapon parameters, 
t 
represents the time 
o .. 
t at which the time dependence f(t) 
is a maximum and N is a 


normalization parameter chosen such that 


00 
(11-2) 
f f(t) dt = 1. 
o 


A graph of this particular time dependence is given in 


Figure 2 using-typical values of a and 
B. 


Transport of Radiation Through the Environment 


After one has an understanding of the nature 
(energy 


and time dependence) of the radiation emitted from a nuclear 


weapon, it is necessary to gain an understanding of how this 


radiation is transported through the surrounding environment. 


That is, if there is a given amount of radiation produced at 


the weapon burst, 
how much of this radiation will reach some 


point which is a 
known distance away. 


If the gamma ray pulse emitted from a nuclear deto- 


nation can be described with respect to time by a function 


f(t) 
as mentioned in the previous section then at any time 


t 
the number of gammas produced per second is given by 
(next 


page) 
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FfGQ~' EMPIRICAL 
TIME DEPENDENCE 
.fOR 
GAMWoA YIELO(THIS 
GRAPH 
IS 
NOT 
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(11-3') 
= 
Y 
f (t) 
Y 
EY 


where Y 
is the gamma yield of the explosion (in MeV), 
E 
is 
y 
y 
the average photon energy 
(in MeV - 
assuming a monoenergetic 


source), and Ny(t) is the' number of gamma rays. 
It can be 


shown 
(see Appendix A) that if the weapon burst is assumed 


to be a point source of gamma's then at a distance R from 


this source the number of unreacted or virgin gammas passing 


through a unit area per unit length of time 
(the virgin 


intensity - 
¢) is given by2 


(11-4) 
¢(R,t) = 


R 
dN 
(t--) 
Y 
c 
dt 


R 
- 6 11 (r)dr 
e 


where 
11 (r) is the total macroscopic cross section 
(see 


Appendix A) 
for photon reaction with the environment and c 


is the speed of light. 
As is evidenced by (11-4), when con- 


sidering high altitude EMP one is concerned only with those 


gammas that arrive at a point R without having been scattered 


or acted upon by the environment in any manner 
(virgin inten- 


sity - 
¢). 
The reason for this is that the high frequency 


EMP associated with high altitude bursts is essentially a 


"line-of-sight" phenomenon. 
That is, if an observer is below 


the burst, then only those electrons that are traveling 


directly on the "line-of-sight" between the observer and the 


weapon burst will contribute significantly to the EMP that 


the observer measures 
(see Figure 3). 
This will be more 


2The R in the argument ofoN 
accounts for time delay 
c 
y 
due to transport of y - pulse over distance R. 
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Figure 3. 
EMP as Line-of-Sight Phenomenon. 
The fields from 
the electron created at 
(B) will contribute to the 
EMP measured at 
(C) whereas the electron created 
-at 
(D) will make no contribution. 


fully explained in Chapter 
11~. 
Although ·it is possible that 


a gamma photon which was not originally on the "line-of-sight" 


could be scattered ohto the "line-of-sight" thereby creating 


an electron which could possibly contribute to the EMP, the 


chances for this are small enough that it will be ignored' 


in the rest of this development. 


Although 
(11-4) 
provides a formalistic solution to the 


problem of the transport of virgin gamma photons, it is de- 


sirable to have a form of 
~(R,t) with which one can obtain 


actual numerical values. 
The main problem with this in 
(11-4) 


is the following integral in the exponential term: 


(11-5) 
R6 ~(r) dr 
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In order to obtain a closed solution bf this integral, it 


is necessary to make .some approximations concerning the 


nature of the function 
~(r). 
First it shall be assumed 


that the macroscopis cross section 
~(r) varies with r only 


as the density of the media through which the gamma photons 


are being transported varies with r. 
In air, 
a good approx- 


imation is to let the density vary exponentially with the 


altitude 
Z above sea level 
(see Figure 4) 
such that 


(11-6) 
p(Z) 
-Z/Z 
= p e 
0 
o 
• 


where p{Z) is the density of air at an altitude Z, 
p 
is 
o 
the density of air at sea level, and Z 
is a scaling factor 
o 
(if Z is measured in kilometers then 
Zo~7). 
If one lets 
~o 


represent the macroscopic cross section at sea level and 


makes the substitution Z = HOB-r then using the approximations 


outlined above 
~(r) becomes 


(11-7) 
11(r) = 
-(HOB-r)/Z 


~ 
~oe 
0 


The integral in 
(11-5) 
can now be evaluated using 
(11-7) 
such 


that 


(11-8) 
R 
f 
~(r)dr 


o 


R 
-(HOB-r)/Z 
= f 
~ e 
0 
o 
0 
dr 


Upon integrating, the right hand side of the above becomes 


-(HOB-r)/Zo 
dr 
-HOB/ZO 
r/Z R 
= 
~ Z e 
(e 
I 
o 
0 
o 
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Figure 4. 
Geometry of High Altitude EMP. 


or 


(11-9 ) 
R 
J 
11 (r)dr 
o 


-HOB/Z 
R!Z 
= 11 
Z e· 
0 
(e 
0 
- 
1) 
o 
0 
(3) 


Substituting 
(11-9) 
into 
(11-4) 
one has a closed solution 


for the virgin intensity of gamma photons at a distance R 


such that 


(11-10) 
<I> (R,t) 
-HOB/Z 
R/Z 


exp{l1 Z e 
o(l-e 
o)} 
o 
0 


Using the expression for dNy given in 
(11-3) 
the virgin 
dt 


3The approximations made in solving this integral were 
obtained through private correspondence with Capt. Terry 
c. 
Chapman, 
GNE-74M, AFIT. 
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intensity can be written as 


(II-II) 


R/Z 
(l-e 
o)} 


The reader should be aware that it has been tacitly assumed 


that the macroscopic cross section 
~(r) is independent of 


time. 
For early time periods following a weapon burst this 


is a safe assumption. 
The air chemistry involved in late 


time calculations concerning 
~(r) are beyond the scope of 


this thesis. 


Creation of Electrons 


There are several ways in which gamma photons can inter- 


act with the molecules of air to produce electrons. 
These 


processes include the photoelectric effect, pair production, 


and the Compton effect. 
The probability of anyone of these 


events occurring to a given gamma photon depends upon the 


energy of that photon and the material through which it is 


traveling. 
This fact 
is illustrated in Figure 5 on the 


following page which is a graph of the mass attenuation 


ff " 
4 
d' 
h 
. 
c('e 
lClents 
correspon lng to t 
e varlOUS processes as a 


function of the energy of the photon in air. 
The curve 


corresponding to total attenuation is the overall sum of 


all the individual effects. 
With respect to Figure 5 one 


can see that for relatively low energies «0.1 MeV) 
the 


photoelectric effect is the predominant factor in the total 


4 h 
'ff' 
. 
T e mass .attenuatl0n coe 
lClent 
cross section 
(see Appendix A) divided 
the photon environment. 
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i s 
(A>~) =(Citj,) 
T (""/1" ) + ( ~~ ) 


where 
0;: 
,7 
and J< 
are the corresponding linear 
coefficients for 
Compton absorption, photoelectric 
absorption, 
and pair production;;::> 
is the air 
density. 
For 
exar;;ple a! 1.5 
~!eV the Compton 
absorption is 
2.4 x 10-'> meter 2/kg of 'alr; 
at 
0 0 
C and 760 
mm Hg pressure, 
the density of air 
is p 
:: 1.293 kg/m 3 . 
(Eef~ 14, p.40) 
18. \ 


GNE/PH/74-7 


attenuation whereas for relatively high photon energies 


(>10 MeV) pairproduqtion provides the major contribution. 


However, for the intermediate photon energies 
(between 0.1 


and 10 MeV) 
the Compton effect is the predominant inter- 


action for photon attenuation. 
Therefore, since weapon 


gamma radiation consists of photons with energies between 


1 MeV and 5 MeV the Compton process is the effect which 


shall be considered with respect to the creation of electrons. 


The Compton effect (scatter) is the process whereby a 


• 
gamma photon is scattered by an atomic electron such that 


the energy of the gamma photon is decreased and the electron, 


which will be assumed to be initially at rest, is ionized 


and set into motion. 
For gamma photons of 1.5 MeV the 


atomic electrons, because of their relatively low binding 


energy, may be· considered as free electrons which are at 


rest. 
The geometry of the situation is illustrated in 


Figure 6. 
If E 
is the incident photon energy then one· 
y 


can easily show (see Appendix A) that the energy of the 


scatter photon is given by the relation 


(II-12) 


I 
E 
= 
Y 


EY 


1.+ 
(1 - cos 


where e is the angle between the direction of the incident 


and scattered photon 
(see Figure 6), mo is the electron 


rest mass, 
and c is the speed of light in a vacuum. 
The 


kinetic energy of the scattered electron T 
can be easily 
e 
found using the conservation of energy 
(assuming elastic 
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Figure 6. 
Geometry of Compton Collision. 
(h is Planck's 
constant and v is the photons' 
frequency) 


collision). 
That is 


(II-13) 
Te = E -E' 
Y 
y 
.(l-cos 8)E jm c 2 
= E { 
yo} 
y l+(l-cos 8)E jm c 2 
y 
0 


As was mentioned earlier, 'the only electrons which are' of 


interest with respect to the high altitude EMP are those 


which are moving on the "line-of-sight" that is, those 


electrons corresponding to ¢ = O. 
This is not entirely 


tYue. 
It will be shown that electrons which vary slightly 


off the "line-of-sight" do contribute to the EMP. 
Although 


the 
ma~ority of electrons scattered by 1.5 MeV gamma photons 


are in the "line-of-sight" direction 
(¢=O) ,one can show 


that there is a distribution of directions for the scat- 


tered e18cb"ons ranging from ¢=Oo to ¢=90o 
(Ref 17, pp 16-9). 


It is possible to account for those electrons scattered off 
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the 
II line··of-sight II but which still contribute to the EMP. 


However, this is a secondary effect and will be ignored in 


this"thesis. 
Therefore it will be assumed that all the 


electrons created by the gamma photons travel in the initial 


direction of the gamma photon 
(~=OO). 
To find Te for these 


electrons, one uses 
(II-13) with 8=180 0 
since momentum must 


be conserved. 
Thus 


(11-14 ) 
Te-Los = 


E 


Y 


where T 
L 
is the kinetic energy for the "line-of-sight" 
e- os 
electrons. 


Since it will be necessary later to superimpose the 


individual effects of these Compton electrons, it is de- 


sirable to know how many of these are created at various 


distances from the burst. 
A quantity which is needed in 


determining this is the reaction rate or the number of 


compton events which. occur in a unit volume over a unit 


time period. 
.If A is the macroscopic cross section for 


the Compton interaction and 
~(r,t) is the virgin intensity 


of gamma photons at a distance r 
from the burst, then the 


reaction rate for Compton scatter of virgin photons is 


given by the relation 


(11-15) 
R (r,t) = A{r) 
~(r,t) 
c 


where Rc is the reaction rate. 
If one assumes that each 


one of these events creates an energetic electron which 


contributes to the EMP {it's motion is on the "line-of- 


sight") 
then the time rate of change of the number density 
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of these "primary" electrons would 
b~ given by the relation 


or 


(11-16) 


dN 
. (r, t) 
prl 
dt 


dN 
. (r, t) 
prl 
dt 


= 
R (r,t) 
c 


= A(r) 
<I>(r,t) 


where N 
. is the primary electron number density. 
Therefore, 
prl 
at any time t after the detonation and at any distance R from 


the point of detonation, the number 
den~ity of primary elec- 


trons is 


t 
N 
. (R, t) 
= A(R) 
f 
<I> (R, t 
I 
) dt I 
prl 
° 


or using equation 
(11-10) 
and 
(11-4) 
-HOB/Z 
R/Z 


exp{~ Z e 
o(l-e 
o)} 
· 0 0 
(11-17) 
N 
. (R,t)=A(R)------------ 
prl 
4TIR 2 


t 
f f(t' 
o 


~ ) dt I} 


To see physically what is represented by (II-I?) 
one 


should understand how the gamma energy is deposited through- 


out the atmosphere with respect to distance from the burst. 


This can be done by considering 
(11-17) 
for late times 


(t+oo ). 
If one assumes a normalized time function then in 


the limit as t 
approaches 


(11-18) 


00 the expression for N 
. becomes 
prl 
-HOB/Z 
R/Z 
Y 
exp{~ 
Z e 
o(l-e 
o)} 
N 
. (R,t+oo ) 
== A(R)EY 
0_0 
_ 
prl 
y 
4TIR 2 


The function 
A(R) 
can be approximated using the exponential 


approximation illustrated in 
equ~tion (11-7). 
A graph of 
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N 
. (R,t+oo ) is given in Figure 7 on the following page where 
prl' 


V , which is the total attenuation cross section, has been 
o 
taken to be equal to the cross section for the Compton inter- 


action. 5 
It should be noted that it is a rather common 


practice in much of the technical literature concerning the 


high altitude EMP to use A for 
~ 
in the above expression 
o 
since the Compton effect is a predominant interaction in 


the energy ranges of interest. 
This is seen to be a valid 


assumption if one refers to Figure 5 which shows that the 


total attenuation coefficient is due almost entirely to 


Compton attenuation for 1.5 MeV gamma photons. 
With regard 


to Figure 7, it is obvious that most of the primary electrons 


are created in a layer of air between thirty and fifty kilo- 


meters. 
For this reason, 
when speaking of the high altitude 


EMP, one often refers to the "absorption layer" which is 


that portion of the atmosphere in which most of the gamma 


energy is absorbed and most of the primary electrons are 


created as described above. 


Summary 


There have been several assumptions and points made 


in this chapter which should be reiterated at this time. 


First, although the gamma photons which emanate from a 


nuclear weapon burst have an energy spectrum ranging from 


around 1.0 MeV to 5.0 MeV, it will be assumed for the re- 


mainder of this study that the gamma spectrum is mono- 


5These values were obtained from the American 
Institute of Physics Handbook 
(Ref. 
5, 
p 8-93). 
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energetic with an energy of 1.5 MeV. 
With respect to the 


creation of primary electrons it has been asserted that 


only the virgin or unscattered gamma photons are important 


for the high frequency EMP. 
Finally, the attenuation of 


the virgin photons with energy of 1.5 MeV has been assumed 


to be due to the Compton effect and all other processes 


have been ignored. 
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III. 
Electromagnetic Theory of the High Altitude EMP 


The purpose of this chapter is to describe and to develop 


the microscopic model. 
Difficult and unwieldy derivations ap- 


pear in the Appendices 
an~ are referenced in the body of the 


chapter. 


The first section of this chapter is devoted to the gen- 


era1 field equations as developed from Maxwell's equations for 


a single electron. 
This is followed by a section concerning 


the equations of motion of an electron in the upper atmo5phere. 


The results of the first two sections are then combined in the 


third section to obtain-the field equations for an electron 


moving in a specific manner. 
The effects of all those elec- 


trons which contribute to the high altitude EMP are then 


summed up in the fourth section. 
The effects of air conduc- 


tivity are ignored in this chapter but are discussed in 


Chapter v. 


General Field Equations for 
~ Single Electron 


The most basic equations in electromagnetic theory are 


Maxwell's equations. 
They are presented here without justi- 


fication or derivation as follows 6 


'i/ x E(r,t) = 
aB (r,t) 
at 


'i/ x B (r, t) = 11 J(r,t) 
+ 11vE:v 
aE(r,t) 


(III-I) 
v 
at 


p (r,t) 


'i/ . E (r, t) 
c 
= 
E:v 


'i/ . B (r, t) = 0 


6The interested reader is referred to the text Electro- 
magnetic Fields and Waves by Lorrain and Corson. 
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where E is the electric field, 
B is the magnetic field, 
p 
is 
c 
the charge density, 
J is the current density, 
~ 
is perme- 
v 
ability of vacuum, 
£ 
is the permittivity of vacuum, 
r is 
v 
the position at which the fields are measured, 
and t is the 


time at which they are measured. 


If one defines the quantities A(r,t) 
and V(r,t) 
such 


that7 


B(r,t) = v x A(r,t) 
(111-2 ) 
E(r,t) = -VV(r,t) 
_ aA(r,t) 
at 


and then manipulates the equations (III-I) 
using the above 


substitutions, the following relations are obtained 
(see 


Appendix B): 


(III-3 ) 
p (r, t) 
c 


(111-4 ) 
2- - 
V2A(r,t) 
- 
~ 
E 
a A(r,t) = - 
~ 
E 
J(r,t) 
v 
v 
at 2 
v v 


- 
- 
aV(r,t) 
+ V{VoA(r,t)} + 
V{~VEV 
at 
} 


Since a vector field is not uniquely defined unless one 


specifies both its curl and divergence and since only the 


curl of A has been specified in the above, one can "choose" 


the following relationship concerning the divergence of A 


(III-5) 
V-A(r,t) = - 
~ 
E 
aV(r,t) 
v v 
at 


This procedure is known as picking the guage"of A and the 


particular 
gaug~ chosen above is called the Lorentz gauge. 


Upon substituting 
(111-5) 
into 
(111-3) 
and 
(111-4) 
one 


7A is often called the vector potential and V the scalar 
potential. 
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obtains the following wave equations8 


p 
(r, t) 
c 
= - 
E:v 


and 


- 
11 
J(r,t) 
v 


which can be written in the following abbreviated manner. 


(\7 2 
- 
1 
3 2 
IA(r,t) I = 


I-~V J'(r ,t) 
I 
(111-6 ) 
-) 
V (r, t) 
- 
1 
p (r, t) 
c 2 3t 2 
- 
C 
E:v 
or 


'" 
IA(r,t)I 
,-11v 
J(r~t) I 
= 
(111-7) 
L 
V (r, t) 
-1 
pc(r,t) 
t,r 


E:v 


The sOlution to equation 
(111-7) 
can be written formally as 


follows: 


(III-8 ) 
IA(r,t) I = Jd 3r'Jdt' G(r,tlr;t') 
V (r, t) 


-11 
J(r;t') 
v 
-1 
pc(r;t') 


E:v 


where 
G(r,tlr~t') is the Green's function associated with the 


operator Lt 
and which must satisfy the following equation 
,r 


'" 
(III-9) 


where 0 is the Dirac delta function. 


In order to obtain a closed form solution to 
(111-8) 
one 


must find the form of the Green's function which is determined 


by (111-9). 
This can be easily done 
(see Appendix B) 
using 


Laplace Transfobn techniques to obtain the following solution 


(III-10) 
1 
oCt' 
= - 
4''!f 


_ {t 
- Jr - r' L}) 
c 


8The relation 
11 
E: 


V 
V = 1 
has been used. 
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Thus, since this source is localized rather than distributed 


through space the current density, as well as the charge 


density, can be represented using the Dirac delta function 


such that 


(III-II) 


and 


(111-12) 


where 13 
v 
= c' 
e is the amount of charge on an electron, and 


x(t) is the location of the electron at time t. 
Figure 8 


illustrates the geometry of the situation. 
Using the above 


expressions for 
J 
and 
p 
and the solution of G(r,tlr;t') 
c 
given in 
(111-10), the equations 
(111-8) 
for A and V become 


(111-13) 


and 


(111-14) 


o(t' -{t- Ir-r' I}) 
c 
\r-r'l 


o'(t'_{t_1r-r' I}) 
c 


The equations given above could be solved to obtain a closed 


form for A and V. 
Then the operations indicated in equations 


(111-2) 
could be performed to find expressions for E and B. 


However, it is less tedious if one simply performs the oper- 


ations mentioned directly on the integral forms of A and V 


given in 
(111-13) 
and 
(111-14). 
Performing these operations, 
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Figure 8. 
Geometry for a Single Electron Source. 
One should 
note that J(r,t) 
and 
p (r,t) are zero unless rtt)=x(t) 
c 


one obtains, after much manipulation 
(see Appendix B), the 


following relations for E(r,t) 
and B(r,t): 


(III-IS) 
E (r, t) 


• 
+ 1 I 
c 


" 
{(n-S) 
x S} Iret} 


and 


(11.1-16) 
B (r, t) 
1 
" 
= -{n xE(r,t)} 
c 


where R=!RI=lr-r' I is the distance between the source 
(electron) 


A 
R 
and the observer, 
n=R is a unit vector indicating the direction 


between the source and the observer, B is the time derivative 


of the "velocity" of the electron (acceleration) 
~ 
and 
K~l-n·~. 
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The bracket with the subscript "ret" relates to the fact that 


the quantity in the brackets is to be evaluated at the re- 


tarded time, t'=t-{R(t')/c}. 
That is, what the observer 


"sees" at time t is that which was created at an earlier time 


t' 
where the difference is the time that is required for the 


signal to traverse the distance R. 
One should note that the 


fields indicated by 
(111-15) 
and 
(111-16) are divided into 


"velocity fields" which are independent of acceleration and 


"acceleration fields" which are dependent on ~ in a 
line~r 


fashion. 
One can show that the velocity fields are essen- 
-2 
tially static fields falling off as R 
, whereas the 


acceleration fields are typically radiation fields falling 


off as R-l 
(Ref 10, p 467). 


For long range effects of the EMP one is concerned 


with the radiated fields rather than the static fields 


since the magnitude of the radiated field represents a flow 


of energy. 
Therefore, the fields of interest for the high 


altitude EMP depend on the acceleration of the Compton 


electrons as well as their velocity. 
That is, in order to 


solve for the radiated fields one needs analytical expressions 


for 
~ (t) . 


Electron Kinematics in the Atmosphere 


Recall that the primary source of the high frequency, 


high altitude EMP is the Compton electrons created by the 


gamma radiation from the weapon. 
Since the radiated fields 


depend on the acceleration of the electron, we shall consider 


two kinds of forces which can act upon the electrons 
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which are created in the absorption layer 
(30-50 Km). 


The first of these is the force which results from the inter- 


action of the electron 
(a charged particle) with the earth's 


magnetic field. 
The result of this interaction is a com- 


ponent of acceleration which is perpendicular to the 


instantaneous velocity of the electron. 
Secondly, the 


electron as it travels through the.atmosphere, undergoes 


many collisions with the surrounding atoms causing it to lose 


energy and consequently to de-accelerated. 
The result of 


• 
these collisions is a component of acceleration which is 


9 
parallel 
to the instantaneous velocity. 
Both of these 


forces will be explained in greater detail in the following 


paragraphs. 


Perpendicular Acceleration. 
The basic equation of 


motion for a charged particle under the influence of an 


electric or magnetic field is-the Lorentz force equation or 


(111-17) 
p = q(E + v 
x 
B) 
. 
where p is the time derivative of the particle's momentum, 


q is the particle's charge, v is its velocity, and E and B 


represent the electric and magnetic field respectively. 
For 


the case of the electron in a magnetic field only, 
(111-17) 


becomes 


(111-18) 
p = - e(v x B) 
(see Figure 9A) 


Using relativistic mechanicslO one can easily prove the 


following identity between the momentum and the total 


9For a zero order approximation 
10One should be aware of the fact that the electrons 
which produce the EMP are energetic enough to be considered 
relativistic 
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Figure 9. 
Electron Moving in a Magnetic Field. 
The 
(B) part 
of the above is the same as the 
(A) part except a 
definite set of coordinate axis has been defined. 


energy 
(~l) of a particle in a uniform magnetic field 
(Ref 15, 


= - e(v x 
B) 


z-axis of a coordinate systen in the 


Using trois relation in 
(111-18) 
one finds the following: 


(111-19) 
~l dv 
-; dt 
·c 
If one defines the 


direction of the B field then 
(111-19) 
can be written in 


terms of the 


(111-20) 


components of the 


- ec 2 B 
V 
=---V 
x 
~1 
y 
ec 2B 
V 
=---V 
Y 
~l 
X 


vectors as: 


V 
= 0 
z 
= Voz 
(constant) 
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/ 
In 
o~der to solve the two coupled equations in 
(111-20) 
one 
_ec 2 B 
can let w=~ and then multiply the first equation by i=;::r- 


and add it to the second 
(Ref 15, 
p 59). 
This yields 


ddt 
(V 
+ iV ) 
= - 
iw(V 
+ iV ) 
y 
~ 
y 
x 


The solution to this is of the form: 


where g is a complex constant. 


(111-21) 
V 
+ iV 
= ge -iwt 
y 
x 
-iao 
One can write g as g=V e 
o 
where Vo and ao are real constants. 
Then using this in 


(111-20) 
and separating real and imaginary parts one finds 


that 


(111-22) 
V 
= -V 
sin(wt), 
x 
0 
V 
= Vo cos(wt) 
y 
where ao which represents the initial phase 
(t=O) has been 


taken as 0 for simplicity. 
From 
(111-21), it is obvious that 


the following relation holds: 


(111-23) 
V 
= Iv 
2 + V 
2 
'= 
Ivl 
sin 
T 
(see Figure 9) 
o 
x 
y 


dx 
Integrating 
(111-21) 
once more 
(Vx=dt) 
one finds that 


(111-24) 
x (t) = x 
+ rlcos(wt), 
y (t) 
= Yo + rlsin(wt) 
0 
where 


V 
V 
~l 


(111-25) 
0 
0 
r l = - 
= 
w 
ec 2B 


Also if one integrates the equation for V2 in 
(111-20) 
the 


result if 


(111-26) 
z(t) = z 
+ V 
t 
. 
0 
oz 


Thus an electron moving in a uniform magnetic field with no 


other external forces will move in a helix 
(see Figure 10) 


having its axis along the direction of the magnetic field and 


with a radius r l 
given by (III-25). 
In the special case of 
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Figure 10. 
Trajectory of an Electron in a Uniform Magnetic 
Field. 
Part 
(B) is an "end on", cross-sectional 
view of 
(A). 


the electron having no initial component of velocity in the 
. 
. 
direction of the B field 
(V 
=0), it moves along a circle in 
02 
the plane perpendicular to the field 
(Ref 15, p 60). 
~he 


acceleration for the electron under these circumstances is 


also constant and can be found from 
(111-20) 
using 
(111-23) 


and 
(111-25). 
That is 
.. 
. 
IvI 2 ="V 
2 + V 
2 + V 
2 = w2 (V 
2 + V 
2) 
= w2V 
2 
X 
Y 
2 
X 
Y 
0 
or 


(111-27 ) 
V 
4 
o 
r 
2 
1 
= illJ sin 
T) 4 
r 
2 
1 


V 
4 
o 


One can see that for an electron of a given energy the 


acceleration will be maximum when the incident electron is 


perpendicular to the B field 
(T=900). 
Thus the sin T term 


is seen to represent a scaling factor, 
the magnitude of 
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which depends on the angle of incidence of the electron 


with the magnetic field. 


Parallel Acceleration. 
When speaking of the "slowing 


down" or "stopping" of an electron in a material such as air 
it is useful to define two quantities known as the stopping 


power and the range. 
The stopping power is defined as the 


amount of energy lost by an electron per unit length of 


path through the stopping material 
(Ref 6, 
p 185). 
The 


stopping power is symbolized by 
-d~/dx since it signifies 


the rate of energy change with respect to distance traveled. 


The range of· the electron, or the length of the path it 


travels before corning to rest, can be defined in terms of 


the stopping power a.s follows: 


(111-28) 
o 
= 
- 
J dx = 


~1 


o 
d~ 


- 
J 
d~/dx 
~1 


where Rl is the range of an electron with initial energy of 


~1. 
If one considers the case where an electron loses an 


increment of energy 
~~ so that its energy becomes 
~2=~1-~~' 


then one can compute the distance 
(Rl ,2) that the electron 


must travel to lose this energy increment with the following 


relation: 


~2 
= - 
J 
d~ 
-:-:::-...,-;:-- 
~1 
d~/dx 


or separating the integration into two intervals 


or using 
(111-28) 


o 
= ( - J 
_d~~_ 


~1 
d~/dx 
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o 
) 
- 
( - 
J 
d~ 
) 


~2 
d~/dx 


(:tII-,29) 
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Rl ,2 = Rl 
- 
R2 


Equation (III-29) is necessary in this development 


because there are no simple, analytical expressions for the 


forces involved in the slowing down of an electron by a 


material. 
In fact, the majority of reliable information con- 


cerning this subject has been derived from experiment and is 


easily obtainable in the form of the total range (i.e. 
Rl or 


R2 and not Rl ,2) of electrons with known initial energy. 


Therefore, in order to obtain information concerning the. 


acceleration of an electron in such a situation it will be 


necessary to construct a numerical model based on physical 


assumptions. 


The first assumption which shall be made in this model 


is that the electron remains on its original path while it 


is being slowed down and is not scattered in other directions. 


That is, it is assumed that the forces acting on the electron' 


are viscous in nature, changing only the energy of the elec- 


tron and not affecting its direction. 
This approximation is 


best for times soon after the creation of the Compton elec- 


tron and loses its validity after the electron has lost a 


large portion of its energy. 
As was mentioned, the forces 


acting upon the electron are assumed to be viscous in nature 


also implying that the slowing down procedure is continuous 


and that the acqeleration function is a smooth curve. 
An 


empirical relation for the range of an electron as a function 


of its energy has been found using these assumptions. 
This 


relation is 
(Ref 17, 
p 
l6-4)~ . 
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2 
= 0.61 ';1 
';1 +0.26 


with 
';1 in .MeV and Rl 
in gm/cm 2 • ll 


In order to numerically approximate the continuous 


slowing down of the electron it will be assumed that the 


electron loses small increments of energy in a point-wise 


manner along its path. 
For example, if the electron has an 


initial energy of ';1' 
and the energy loss increment is 6';, 


then it will be assumed that the electron travels a distance 


Rl ,2' given by 
(I11-29), with a constant energy ';1. 
After 
it has traveled this distance the electron then instanta- 


neously "loses" 
6'; amount of energy thereby acquiring an 


energy of ';2=';1-6';. 
This procedure is then repeated in a 


step-wise fashion throughout the slowing of the electron. 


Figure 11 gives a pictorial explanation of the situation. 


It should be pointed out that the accuracy of this approach 


increases as the magnitude of 6'; decreases. 
One should also 


note that the successive distances between energy increments 


(i.e. 
Rl ,2' 
R2 ,3' etc) 
are not necessarily equal. 


The acceleration can be approximated by giving the 


acceleration at the ith point 
(see Figure 11) 
such that 


(III-3l) 
{Vi +l 
- Vi_I} 
{Xi +l - Xi-I} 
v. = 
1. 
x i +l - x. 
1 
t i +l - t. 
1 
1.- 
1.- 


where V. is the velocity corresponding to .;. which the 
1. 
1. 


llwhen measured in these units Rl 
is simply the actual 
distance traveled in cm multiplied by the density of the 
stopping material. 
This is convenient for experimental 
purposes. 
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Figure 11. 
Approximation of Continuous Slowing Down of Electron. 
The i th point is not at the point of energy change in 
order to avoid the question of infinite acceleration. 


particle has at the i th point, x. is the total distance from 


~ 
the point i=o to the i th point, 
and t. is the time required 
. 
~ 


to traverse the distance 
x .• 
If one realizes that the: second 


~ 


term on the right hand side of 
(111-31) 
approximates the ve- 


locity which the electron has at the i th point, then the 


expression for V. 
can be written as: 


.~ 


(111-32) 


Therefore, 
one can use 
(111-32) 
for an electron of initial 


energy 
~ by incrementing the energy as indicated above and 


using 
(111-29) 
and 
(111-30) 
for the values of x. while sub- 


~ 


stituting the values for V. corresponding to each energy' 
~ .• 


~ 
~ 


This will yield a series of points for V from which a continuous 
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function can be implied. 


In order to illustrate the method described above, the 


calculations were performed for 1 MeV electron at an altitude 


of 40 Km. 
Figure 12 illustrates the form of the acceleration 


by showing the relative velocity 
(S=V/c) 
as a function of 


time. 
Although the actual numbers are different for various 


altitudes, the general shape of all the curves is the same 


for the different altitudes. 
One should note that the major 


portion of the electron's acceleration occurs in a relatiyely 
• 


short time at the end of the electron's trajectory. 


Electron trajectory. 
The electron in the atmosphere under- 


goes both the parallel and perpendicular acceleration simulta- 


neously. 
Thus, the developments of the previous two sections 


need to be combined in some manner in order to give the elec- 


tronls true trajectory. 
Since it is not possible to derive a 


self-consistent simple force equation which would simulta- 


neously consider both effects, it is necessary to again create 


a numerical model based on physical assumptions. 


It will be assumed that the effect of the air molecules 


(parallel acceleration) 
on the electron is the same as given 


in the previous section. 
That is, the energy of the electron 


changes by a definite increment at various points along its 


path. 
In between these changes, the electron moves at a con- 


stant energy and consequently at a constant speed. 
However, 


since there is a uniform magnetic field 
(perpendicular accel- 


eration) acting upon the electron, the electron does not 


traverse a straight line in between energy changes. 
Rather, 


the path forms a segment of a helix, the radius of curvature 
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I 
I 
of 
wh~ch is given by equation (III-25) 
since the speed of 


the electron 
(and energy) is constant over this path. 
The 


arc length or actual distance the electron travels for each 


segment between energy changes is given by equations (III-29) 


and (III-30). 
As will be shown later, an important variable 


is the angle which the instantaneous velocity vector of the 


electron makes with the initial direction (line-of-sight) of 


the electron. 
If one measures this angle at the i th point 


(as defined in the previous section) 
then it can be 
give~ in 


terms of previously defined quantities as follows 
(see 


Figure 13): 


L 
R. 
1 
. 
i=l 
1- ,1 
(111-33) 
l/J. 
= 
1 
L 
r. 
1 
i=l 
1- 


h 
,I.' 
th 
1 
d 
t 
h 
.th 
. 
were 
0/' 
1S 
e ang e as measure ate 1 
p01nt. 
1 
To give 


the reader a picture of the electron's trajectory using this 


model, calculations were made for the various quantities 
12 
defined at 40, 
30, and 20 Km. 
These results are summarized 


in Table 2. 
One should note the differences in ranges for 


the various altitudes. 
It should also be pointed out that the 


electron tends to spiral inward and at higher altitudes will 


make several complete revolutions before coming to rest 
(see 


l/J. 
column in Table 2). 
1 


varies exponentially (i.e. 
p=p e 
o 


12These calculations assumed that the 
-z/z 
0) • 
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Table 1. 
Calculations for Electron Trajectory. 
(I13I=10- It weber/m 2 ) 
Electron 
Total Path 
Total Time 
Radius of 
1/1. 
Energy 
Traveled 
Elapsed 
Curvature 
1. 


(MeV) 
(meters) 
(shakes) * 
(meter) 
(radians) 


ABSORPTION ALTITUDE = 40 
KM 


1.00 
0.0 
0.0 
47.49 
0.0 
0.95 
73.94 
26.19 
45.71 
1.56 
0.90 
147.88 
52.51 
43.92 
3.17 
0.85 
221.83 
78.95 
42.13 
4.86 
0.80 
295.77 
105.55 
40.32 
6.61 
0.75 
367.25 
131.42 
38.50 
8.39 
0.70 
441.19 
158.39 
36.67 
10.31 
0.65 
512.67 
184.67 
34.82 
12.26 
0.60 
581.68 
210.29 
32.95 
14.24 
0.55 
653.17 
237.13 
31.05 
16.41 
0.50 
719.72 
262.45 
29.14 
18.55 


ABSORPTION ALTITUDE = 30 
KM 
1.00 
0.0 
0.0 
47.48 
0.0 
0.95 
17.72 
6.28 
45.70 
.37 
0.90 
35.44 
12.58 
43.92 
.76 
0.85 
53.16 
18.92 
42.13 
1.16 
0.80 
70.88 
25.29 
40.32 
1.58 
0.75 
88.01 
31.50 
38.50 
2.01 
0.70 
105.73 
37.96 
36.67 
2.47 
0.65 
122.86 
44.26 
34.82 
2.94 
0.60 
139.40 
50.40 
23.95 
3.41 
0.55 
156.53 
56.83 
31.05 
3.93 
0.50 
172.48 
62.90 
29.14 
4.45 


ABSORPTION ALTITUDE = 20 
KM 
1.00 
0.0 
0.0 
47.48 
0.0 
0.95 
4.25 
1.50 
45.70 
.09 
0.90 
8.49 
3.02 
43.92 
.18 
0.85 
12.74 
4.53 
42.13 
.28 
0.80 
16.99 
6.06 
40.32 
.38 
0.75 
21.09 
7.55 
38.50 
.48 
0.70 
25.34 
9.10 
36.67 
.59 
0.65 
29.44 
10.61 
34.82 
.70 
0.60 
33.41 
12.08 
23.95 
.82 
0.55 
37.51 
13.62 
31.05 
.94 
0.50 
41.34 
15.07 
29.14 
1.07 


*1 shake=10- e second 
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. .Fields for Electron Moving in.SpecificManner 


The general field equations for an accelerated electron 


have been derived and are given by equation 
(111-15). 
Further- 


more, it has been asserted that the electron is acted upon by 


a force which is perpendicular to its instantaneous velocity 


(B 
- 
field) 
and a force which is parallel to this velocity 


(interaction with air molecules). 
Since some insight has been 


gained concerning the nature of the acceleration of the elec- 


tron when subjected to these forces, it is now possible to 


arrive at relations for the fields which will be radiated by 


the electron as it is accelerated in this manner. 


First, however, it is useful to define a quantity known 


as the Poynting vector S which represents the energy which 
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flows across a unit area per unit 
ti~e and is 
associate~ 


with radiated 
fields~ 
If 
E 
and B. represent the radiated 
a 
a 
portions of the fields then 5 is defined as 


(111-34) 
5 = 


or using 


5 = 


(111-16) 


1 
- 
-- {E 
x 
C].lv 
a 


One can easily show that E 
is perpendicular to n which 
a 
implies that 


(111-35 ) 
5 = 


• 


Using 
(111-35) 
one can then find the energy which is flowing 


in the direction of the observer by the following 


5 
0 


or, using the radiation term in (III-IS) 


(111-36) 


A 
• 


{(n-S) x S}1 2 }ret 


Using 
(111-36) 
one could then compute the amount of energy per 


unit area radiated during a finite period of acceleration, t' 
= 


Tl 
and t' = T2 , and detected in a corresponding interval, 


t = Tl 
+ {R(Tl)/C} and t = T2 + {R(T2 )/C}, by the observer 


with the following: 


T2+R(T2 )/C 
A 


U = f 
{5 
o n}ret dt 
Tl+R(Tl)/C 


or 
T2 
A 
(~~ I) 
(111-37) 
u = f 
(5' n) 
dt' 
T1 
45 
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wher~ U is the energy measured. 
It is obvious from 
(111-37) 


that the quantity 
(S'~)(~:I) represents the power 
(energy/time) 


radiat~d per unit area by the electron in terms of the "local 


time" of the electron rather than the retarded time defined 


earlier. 
Thus, with the definition of the solid angle in 


mind, one can define the power radiated per unit solid angle 


and measured by the observer at a distance R from the source 


as follows: 


dP (t' ) 


dn 


or using the relation for retarded time 
(see Appendix B) 


dP(t') 
- 
KR2(S'~) 
dn 
- 


which, using the definition of 
K and equation (111-36), 


becomes 


(111-38) 
dP (t' ) = 
dn 


A 
A 
• 
2 
{In x 
{(n-~) x_B~ ,_ } 


(1 
- 
n 
• 
(3) 


One can reduce 
(111-38) to a simpler form for the case 


corresponding to an acceleration which is acting parallel to 


the instantaneous velocity by recognizing that 


A 
A 
• 
A 
/\_ 
• 
A 
A.. 


(III-39) 
n x {en-B) 
x B} = n x(nxS-SxS)= n x 
(nxS) 
. 
Since S is parallel to S. 
Using 
(111-39) 
and letting 
~ be 
. 
"'- 


the angle between S, Band n 
(see Figure 14 A) 
then one 


finds that 


dPL (t' ) 
• 
2 
sin 2 
(III-39) 
1 
(eS) 
ljJ 
= -- 
~ 5) 
dn 
£ 
c 
4n 
(1-13 cos 
v 
where PL signifies the power radiated which is due only to 


parallel acceleration. 
A similar reduction in form is 


possible for the case where the acceleration is perpendicular 
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to 
t~e instantaneous velocity. 
If one defines n as the angle 


• 
A 


between 6 and the projection of n into the plane perpendicular 


to S (see Figure 14 B) 
then one finds that 


dP (t') 
. 2 
sin 2 
cos 2 
(111-40) 
c 
1 
(e-B) 
1 
{I 
1/1 
n 
} 
= 
dst 
CE: 
4n 
(1-6 cos 
1/1) 3 
y2 (l-6 cos 1/1)2 
v 


where y2=(1_6 2)-1 
and P 
signifies the power radiated which is 
c 
due only to perpendicular acceleration. 


With regard to the high altitude EMP the next step is to 


investigate the relative contributions of radiation due to 


parallel acceleration and perpendicular acceleration. 
Cal- 


culations can be made for 
(111-39) 
and 
(111-40) at various 


points along the electron's trajectory using the model out- 


lined in the last part of the previous section to obtain 


values for the acceleration terms at various points. 
If one 


assumes that the observer is situated on the original line-of- 


sight, then n=180° and 
1/J at each point is the angle 
1/1. defined 
J. 
in the previous section. 
One should note that n=Oo or 180 0 


corresponds to the minimum value for 
(111-40) with respect 


to n. 
The calculations indicated here were performed by the 


author for Compton electrons created at various altitudes 


ranging from 60 to 20 Km. 
The results at 40 
Km for early 


times are presented in Table 2 on Page.49. 
As one can 


readily see, the radiation resulting because of the electron's 


interaction 
wit~ the magnetic field 
(dP /dst) 
far outweighs 
c 
that resulting from the linear slowing down of the electron 


(dPL/dQ). 
Although Table 2 shows only the calculations at 


early times, 
and at one altitude, the results are found 
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Figure 14. 
Geometry for Angular Distribution of Electron 
Radiation. 
Part A is for parallel acceleration; 
Part B concerns perpendicular acceleration. 


to be the same at late times and at all altitudes. 
Thus, 


based on these 
calc~lations, one can say that the electro- 


magnetic radiation which constitutes the high 
altitud~ EMP 


results from the turning of high energy Compton electrons by 


the earth's magnetic field. 


It is beneficial to investigate even further the nature 


of the radiation resulting from the turning of the electron. 


Table 2 shows that for a period over which the radiated power 


drops four orders of magnitude, the velocity, S, of the 
ele~- 


tron is 
essenti~lly constant to 
tw~ significant figures. 
Thus, 


although much time has been devoted to the slowing of the elec- 


. 
h 
h 
13 .. 
h 
h 
. h 
d 
tron 1n t 
e atmosp ere 
, lt 1S seen 
ere tat, Wlt 
regar 


13This development will be very important later, however, 
when discussing air conductivity. 
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Table 2 
dPL 
dP 
Comparison of dn 
and dQc for Electron at 40 Km. 


Electron 
Energy 


(MeV) 
1.0000 
.9964 
.9929 
.9893 
.9857 
.9821 
.9786 
.9750 
.9714 
.9679 
.9643 
.9607 
.9571 
.9536 
.9500 
.9464 
.9429 
.9393 
.9357 
.9321 
.9286 
.9250 
.9214 
.9179 
.9143 
.9071 


f3 
(-) 
.941 
.941 
.940 
.940 
.940 
.940 
.939 
.939 
.939 
.938 
.938 
.938 
.937 
.937 
.937 
.936 
.936 
.936 
.935 
.935 
.935 
.934 
.934 
.934 
.933 
.933 


Total Time 
Elapsed 
(shakes) 
0.00 
1. 82 
3.64 
5.46 
7.28 
9.10 
10.92 
12.75 
14.57 
16.39 
18.21 
20.03 
21.85 
23.67 
25.49 
27.32 
29.14 
30.96 
32.78 
34.60 
36.43 
38.25 
40.07 
41. 89 
45.54 
47.36 


1jJ.1 
(radians) 
0.000 
.108 
.217 
.325 
.434 
.544 


. 653 
.763 
.873 
.984 
1.094 
1.205 
1.316 
1. 428 
1.539 
1.651 
1.764 
1.876 
1.989 
2.102 
2.216 
2.329 
2.443 
2.558 
2.787 
2.902 


dPL 
dn 
(*) 
o 
1.75E-39 
2.24E-39 
1.17E-39 
4.56E-40 
1.69E-40 
6.46E-41 
2.62E-41 
1.14E-41 
5.27E-42 
2.59E-42 
1.34E-42 
7.30E-43 
4.14E-43 
2.43E-43 
1.47E-43 
9.10E-44 
5.77E-44 
3.73E-44 
2.41E-44 
1.58E-44 
1.03E-44 
6.58E-45 
4.10E-45 
1.29E-45 
5.21E-45 


dPc 
dn 
(*) 
1.07E-13 
9.32E-14 
4.81E-14 
7.00E-15 
6.80E-17 
2.98E-16 
. 
4.74E-16 
3.89E-16 
2.69E-16 
1.76E-16 
1.15E-16 
7.64E-17 
3.63E-17 
2.61E-17 
1.93E-17 
1.47E-17 
1.14E-17 
1.10E-18 
7.42E-18 
6.20E-18 
5.28E-18 
4.60E-18 
4.08E-18 
3.70E-18 
3.21E-18 
3.08E-18 


*These calculations were performed using CGS-Gaussian units 
rather than MKS units. 
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to 
th~ high altitude EMP, the radiating electron can be con- 


sidered as moving with a constant speed and being turned by 


the magnetic field of the earth for the time periods of inter- 


est. 
Making this assumption, one can then calculate the angle 


that the instantaneous velocity vector makes with the initial 


direction of the electron 
(the "line-of-sight" direction for 


example) at any time t 
as follows: 


(111-41) 


where Wis the indicated angle and r S is the radius of curva- 


ture given by 
(111-25) corresponding to an electron with 


speed S. 
Using 
(111-41) 
in 
(111-40) 
one could obtain in- 


formation concerning the power or the magnitude of the E 


fields which an observer located along the initial direction 


of the electron would measure as a function of time. 
Such 


a calculation is illustrated in Figure 15 which is a graph 


of normalized values of the time dependence indicated by 


(111-40) 
and 
(111-41). 
An even more intuitive understanding 


of the distribution of this radiation is gained from Figure 16 


which illustrates an iso-field line for an electron with a 


certain instantaneous velocity. 
From these two figures, 
one 


can see that a large portion of the electron's radiated en- 


ergy is directed forward in the direction of its instantaneous 


velocity. 
Also, it is evident that the strength of the fields 


drops off rapidiy as one moves away from this direction. 
Thus, 


an observer measures an initial maximum field strength which 


drops off drastically after the electron turns only a few 


degrees. 
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The Net Effect of Many Electrons 


Line-of-Sight Effect. 
Since there are many Compton 


electrons created in the upper atmosphere by a high altitude 


nuclear detonation, it is necessary to know which of these 


contribute to the high frequency EMP measured by an observer 


at a particular location. 
It has been previously asserted that 


the Compton electrons which contribute to the high frequency 


(early time) 
EMP are those which are created by gamma photons 


moving on the "line-of-sight" direction for the observer.with 


respect to the weapon detonation 
(see Figure 17).14 
There are 


two primary reasons for this which involve the time phasing 


of the radiated fields and the angular distributions of the 


radiated energy. 


The important factor with regard to the time phasing 


aspect is that the Compton electrons are traveling at essen- 


tially the speed of light during the time they radiate the 


major portion of their energy. 
The importance of this fact 


can be understood by considering two gamma photons emitted 


at the same time and moving together on the line-of-sight, 


each of which creates a Compton electron but at different 


times 
(see Figure 17). 
Since the first electron created is 


essentially moving at the same speed as the remaining gamma 


photon, the fields which are radiated in the time interval 


T1 
(Tl 
small enough so that the electron has not appreciably 


turned off the line-of-sight) 
and which propagate along the- 


l4The reader is reminded that it has been assumed that 
all Compton electrons initially move in the incident direction 
of the gamma photon. 
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Figure 16.. Radiation Pattern for Electron in Circular Motion. 
(Note that the coordinate system which defines ¢ 
(and n) is constantly changing since the z-axis 
is defined by S(t).) 
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Time Phasing for-EM Fields for-High altitude EMP. 
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Sp~l, then the two radiated fields indicated 
above will arrive at the observer essentially in 
phase with respect to time. 
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line-of-sight constitute a pulse" the "leading edge" of which 


is almost even with the "leading edge" of the gamma photon. 


Similarily, since the electron created later by the remaining 


gamma photon is also moving with approximately the speed of 


light, the fields which it radiates in the time interval Tl 
and which propogate along the line-of-sight constitute a 


pulse the "leading edge" of which is almost. perfectly matched 


up to the first pulse. 
These two pulses, then, will arrive at 


the observer in such a manner as to result in almost perfect 


constructive interference. 
One can infer from this argument 


that the more the electron's instantaneous velocity vector 


deviates from the line-of-sight the less the component of 


velocity in the line-of-sight direction becomes" and the longer 


the distance which the radiated energy. must travel to the ob- 


server becomes. 
The result of this is that the contribution 


of the radiated pulse to that which the observer measures 


becomes less and less significant as compared to the con- 


tribution of the line-of-sight electrons until it is nothing 


more than background. 
Thus, 
when speaking of "line-of-sight" 


electrons with respect to time phasing one included electrons 


created within a small cone centered about the line-of-sight 


such that their valocity component parallel to the line-of- 


sight is approximately equal to their total velocity and 


such that the distance which the radiated energy travels to 


the observer is approximately equal to that distance which 


the energy radiated by actual line-of-sight electrons travels. 


Independent of the time phasing and more intuitively 
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obvious is the effect of the angular distribution of the 


radiated fields. 
From Figure 15 and Figure 16 one can see 


that'the electron which contributes the most to the power 


the observer measures is the one which is moving directly 


toward the observer. 
Referring to Table 2 one also can see 


that when the angle between the instantaneous velocity and 


the direction to the observer is just a 
few degrees, the 


radiated fields measured by the observer have dropped by 


.. 
an order of magnitude. 
Thus, the electron which is created 
. 
by a gamma whose incident direction is just a 
few degrees 


from the line-of-sight direction will not significantly 


contribute to the early time EMP. 
It can be argued that an 


electron with initial direction away from the observer can 


be turned by the earth's magnetic field such that is eventually 


"points" at the obserVer. 
This argument is valid insofar as 


the motion of the electron is' such that it satisfies the time 


phasing arguments given in the previous paragraph. 


Although it has been shown that the important contribution 


to the high altitude, high frequency EMP is made by line-of- 


sight electrons, one should realize that, when speaking of the 


electrons which contribute to the fields measured by an ob- 


server at a particular location, 
some electrons which are 


created by gamma photons with directions slightly off the 


line-of-sight are to be included. 
The degree of deviation 


from the line-of-sight is determined by the constraints men- 


tioned in the two previous sections. 
It is these qualifications 


which should be kept in mind when speaking of the high altitude, 
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high 'frequency EMP as a "line-of-sight" phenomenon. 


Superimposing Effect. 
When computing the net effect of 


all the electrons which contribute to the EMP one is primarily 


concerned with the total field strength which the observer 


measures as a function of time. 
The variation in time of the 


field strength will depend, 
among other things, on the time 


nature of the gamma photon output. 
From equation 
(11-10) 
in 


Chapter II it has been shown that the number of gamma photons 


passing through a unit area at some distance R from the burst 


and at some time ttl is given by 


(111-42) 
<'P(R,ttl) 
'ty 
= 
f(T")F(R) 
Ey 
R/z 
(l-e 
0) } 
-HOB/zo 
= exp{llozoe 


4nR 2 
F(R) 


where 


and 


Ttl = ttl 
_ Rc 


Furthermore, using 
(11-16) 
one can write the number of Compton 


electrons which are created per unit volume per unit time at 


time til as 


dN 
. (R til) 
prl 
' 
dt 


y 
= EY A(R)F(R)f(T 
tI 
) 
Y 


The number of these electrons which would be created in a 


time interval dt can be found from the following: 


(111-43) 
dN 
. (R til) 
prl 
' 


y 
= EY A(R)F(R)f(T")dt 
Y 


If E (r,t ) represents the general time dependence of the 
e 
e 
field strengths radiated by a single electron and measured by 
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an observer at a distance r 
away and at some time t 
15 
then 
e 
' 
for an electron 
~reated at til the fields measured at time t 


after the weapon detonation would be given by E (r,t-t"). 
e 
Thus an observer would measure the following for 
(V·dN 
.) 
prl. 
electrons created a distance r 
away from him and R from the 


burst: 


E 
(r t) 
= V·dN 
. (R t")E 
(r t-t") 
p 
, 
prl.' 
e' 
.. 


where V is the volume enclosing the electrons which contribute 


as earlier described. 
If the time origin is shifted sucb that 


t=O when T"=O, then the above can be written as 


E (r,t) 
= V{dN 
. (R,T")E 
(r,t-T") 
p 
prl. 
e 


or using 
(III-43) 
and recognizing that the r 
dependence for 


Ee is l/r for radiated fields, 


VY 
E (r,t) 
= --:f.. A(R)F(R){f(T")E 
(t-TI)dT"} 
P 
rE 
e 
y 


EpCr,t) 
just accounts for those electrons created within an 


infinitesimal time increment of dT. 
There will be electrons 


created before and after this increment which will contribute 


to the total field strength, ET(r,t), measured by the ob- 


server. 
Thus, in order to obtain the total field strength 


caused by Compton electrons created within a finite volume16 


centered about a point some distance R from the burst and 


measured at some time t 
by an observer located a distance r 


l5These fields were created at tl<t _£ (retarded time) 
- 
e 
c 


l6For further information concerning this volume see 
Chapter IV. 
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from the electrons, one should integrate 
(111~44) from' 


T=O to T=t-r . 
That is, 
c 


(111-45) 


r 
VY A(R)F(R) t-c 


y 
J 
f(T")E 
(t-T")dT" 
Ere 
y 
0 


As has been indicated, equation 
(111-45) is for the 


electrons created within a certain volume centered about a 


point some distance R from the burst. 
One can sum over all 


such points by letting V=A'dr where A is the cross-sectional 


area of V lying in a plane perpendicular to the line-of-• 
• 
sightl7 and dR is an infinitesimal length parallel to the 


line-of-sight, then integrating 
(111-45) over the absorption 


layer. 
Letting A=A(R) 
one can then write the following for 


the fields which an observer located a distance X from the 


burst would measure at time t. 


(111-46) 
E(X,t) 
Y 
X 
= 
Y J dRA(R)A(R)F(R) 
Ey 
0 
iX - 
R) 


t_(X-R) 
c 
J dT"f(T")E 
(t-T") 
o 
e 


One should note that it was not necessary to make a correction 


for time lags between altitude because of the time phasing 


which was mentioned previously. 


17 h 
. 
t' 
d' 
1 
. 
, 
T e constralnts men lone 
ln tle prevlous sectlon 
regarding the line-of-sight effect determine the size of A. 
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IV. 
Sample Calculations 


The purpose of this chapter is to use the theory devel- 


oped in the previous chapters and obtain data regarding the 


nature of the high altitude, high frequency EMP. 
The tech- 


niques used to solve the various equations and arrive at the 


physical quantities mentioned previously will also be demon- 


strated. 
For simplicity, all these calculations involve an 


observer directly beneath the burst and who is located at 


the magnetic equator. 
Thus the incident velocity of the' 


Compton 
elec~rons will be perpendicular to the earth's 


magnetic field. 
Furthermore, it shall be assumed that the 


earth's magnetic field is uniform over the absorption layer 


-4 
and has a field strength of I.OxIO 
webers/meters. 


Line-of-Sight Area 


As was mentioned in the previous chapter, the high 


frequency EMP is "almost" a line-of-sight phenomenon. 
That 


is, there are electrons just slightly off the line-of-sight 


which contribute significantly to the EMP measured by an 


observer. 
The first task in computing the EMP, then is to 


determine the number of the electrons which contribute to 


the observer as a function of distance from the burst. 
For 


these calculations, it shall be assumed that all the elec- 


trons which contribute, do so in exactly the same manner 


as the electron which is precisely on the line-of-sight. 


Thus, the cross-sectional area located at a certain distance 


R from the burst which contains all the electrons that con- 
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trib~te will be considered to act as a point source to the 


observer. 
With regard to this approximation, it is obvious 


that the effect of 
~lectrons which originate with a velocity 


slightly off the line-of-sight yet which are turned so as to 


"point ll at the observer will be underestimated. 
Similarly, 


the effect of the electrons originating with a velocity off 


the line-of-sight and which are turned away from the observer 


will be overestimated 
(see Figure 18). 
To obtain a measure 


of the error introduced by this, equation 
(111-40) with n=1800 


was integrated bwtween an angle 
~o and 
~o + 6° for 
~o ranging 


between -1° and 1°. 
The per cent difference between each 


integration and the integration between 
~ =0° and 
~ 
+ 6° 
o 
0 
was then computed. 
These results are presented in Table 3. 


For the purpose of this study, it was decided that the error 


introduced by 
~ 
= ± 0.5 was tolerable. 
Thus, the maximum 
o 
angle between the incident direction of the electron and the 


line to the observer in order for the electron to be con- 


sidered as a line-of-sight electron will be set at ±0.5°. 


If one assumes that the cross-section area A(R) at each 


altitude is a circle, then one must find the radius of this 


circle in terms of 
~0(=±0.5°) at each altitude. 
Using the 


geometry indicated in Figure 19, with the observer at sea 


level and letting XLOS be the sought after radius, one finds 


that 


(IV-l) 
XLOS(R) = 
HOB 
+ / 


2tan~o 
( 
HOB 
2 
~~-) 
+ RZ 
2tan1J!o 
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Table 3 


Difference of Total Power 


Received for Electrons 


off Line-of-Sight 


l/Jo 
% 
l/Jo 
% 


DEGREES 
DIFF 
DEGREES 
DIFF 


-1.0 
11. 2 
0.1 .. 
1.3 


-0.9 
10.3 
0.2 
2.7 


-0.8 
9.3 
0.3 
4.0 


-0.7 
8.2 
0.4 
5.4 


-0.6 
7.2 
0.5 
6.9 


-0.5 
6.1 
0.6 
8.3 


-0.4 
4.9 
0.7 
9.4 


-0.3 
3.8 
0.8 
10.6 


-0.2 
. 2.5 
0.9 
11.1 


-0.1 
1.3 
1.0 
12.4 
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Figure 18." 
Iso-Power Lines for Electrons in Three Different 
Configurations Relative to Observer. 
In a set 
time period during which each electron turns the 


.. same distance and consequently sweeps out the same 


number of degrees 
(at least >o/AO) 
then the con- 
tribution of A will be greater than B which will 
be greater than C. 


Therefore, for a circular cross-sectional area A(R) one has, 


. 
{-HOB 
A(R) 
= 4w 
2tano/ 


o 
+ /. ( 
HOB) 2 + 
RZ 
} .2 
2tan% 


A graph was made for XLOS as a function of altitude for a 


weapon detonation occurring at 100 KID and with 0/0=0.5°.. 
This 


is given in Figure 20. 


The assumption that the electrons enclosed in the circular 


area A(R) 
above the observer and centered on the line-of-sight 


act as a point source lying ort the line-of-sight ignores the 


n dependence in equation 
(111-40). 
For the electron orig- 


inating on the line-of-sight, n is always equal to 0° or 180° 
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Figure 19. 
Geamerty for Finding the Line-of-Sight Radius. 
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FIG.20. LINE-OF-SIGHT RADIUS AS 
A FUNCTION OF ALTITUDE. 
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atid cos 2n=1. 
However, there are electrons included in A(R) 


which have values of n ranging from 0° to ±180 0 
(see Figure 


21). 
By saying that the disk with area A(R) 
appears as a 


point to the observer, one is implying that the distance 


XLOS cannot be revolved. 
Therefore, as far as the observer 


is concerned, the electrons all appear at a point and have 


a value of n=O° or 180 0 depending on the orientation of the 


magnetic field. 
To test the validity of this assumption one 


can compare the magnitude of XLOS(R) at various altitudes to 


the magnitude of the distance to the observer. 
For example, 


using Figure' 20 one can see that at an altitude of 50 Km the 


value of XLOS is about 0.22 Km. 
The ratio between this and 


the distance to the observer is .004 so that one could say 


that XLOS«z where z, in this case, is the distance to the 


observer. 
Similarly, at 20 Km the ratio is .007. 
Therefore, 


the assumption of a point source is a relatively good one. 


One should note that the A(R) 
which meets the point source 


specifications falls well within the line-of-sight require- 


ments mentioned in the previous chapter. 


Integration Over Absorption Layer 


The integration over R in equation 
(111-46) 
can be 


approximated by summing over segments of the absorption layer 


each centered about some distance R. 
and having a depth of 


1 


~R. 
The quantity 
~R must be such that the macroscopic cross 


sections and all other quantities depending on R are constant 


over the interval, since all 
q~antities in this interval will 


be assumed to have the value given at R.. 
It was found that 
1. 
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Figure 21. 
Top View of Electron in Line-of-Sight Area A(R) . 


6R=1 Km was a relatively good approximation. 
Using this one 


can now compute the magnitude of the fields resulting from a 


single electron at R. 
which an ovserver would measure at a 
1 


certain time, multiply that number by the Compton electrons 


per unit volume formed at R. at an appropriate earlier time, 
1 


and then multiply by A(R.)6R to get the total effect of "a 


1 


cylindrical volume of electrons centered at R.. 
This pro- 


1 


cedure could then be repeated for all other R. 
in the absorp- 


1 


tion layer and the various contributions summed up in order 


to obtain the net effect. 


Discussion of Results 


The calculations which have been previously discussed 


were performed for a nuclear detonation occurring at 100 
KID 


with a gamma yield of one kiloton. 
The magnitude of the E 
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fields which an observer, located directly beneath the burst 


at sea level and at the magnetic equator, would measure, based 


on these calculations, is presented in Figure 22. 
The length 


-8 
of the pulse is about 5 shakes 
(1 shake = 10 
sec) 
and reaches 


peak fields of the order of 12.0 x 10 4 volts/meter. 
The 


validity of the section of the graph for times after ten 


shakes is doubtful since by this time many of the electrons 


have slowed considerably below the value of B=0.94 which was 


assumed for the calculation. 
In a paper by Karzas and Lptter 


(Ref 11) using a macroscopic approach with similar weapon 


parameters, peak field values of 6.0 x 10 4 volts/meter are 


reported. 
Thus, the microscopic approach does give answers 


which are in the same order of magnitude as accepted theories. 


The difference between the two could be accounted for by the 


fact that conductivity has been ignored here but was accounted 


for in the Karzas and Latter paper. 


Although the previous calculations were done for the case 


where the electron's incident velocity is perpendicular to the 


earth's magnetic field, it is a simple matter to show the ef- 


fects of different orientations using this microscopic model. 


Equations 
(111-35) 
and 
(111-40) 
indicate that the magnitude of 


the E fields measured by an observer are proportional to the 


acceleration, S. 
From equation 
(111-27) 
one can see that for 


an electron of a given energy the acceleration is proportional 


to sin 2 T where 
T is the angle between the incident direction 


and the magnetic field. 
Thus, the magnitude of the fields 


measured by an observer will vary as sin 2 T according to his 
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FIG.22. 
EMP MEASURED 
AT 
GROUND AS 
A FUNCTION Of TIME, 
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magnetic bearing with respect 'to the weapon burst. 
This 


point is illustrated by Figure 23 on the next page. 
The ob- 


server at point E will measure the maximum field magnitudes 


since the angle, 
T, between the initial direction of the 


electrons and the magnetic field is 90°. 
However, the ob- 


server at point A will not measure any fields since the 


direction of the electrons' motion is parallel to the magnetic 


field 
(T=OO). ,The field strengths measured by the observers 


at points B, 
C and D increase as sin 2 T according to the 


respective values of 
T in each case. 


One should be aware of the potential range of the high 


altitude EMP which can be inferred from its "line-of-sight" 


nature and the facts of the previous paragraph. 
'Since the 


high altitude, high frequency EMP is essentially a line-of- 


sight phenomenon, the effects of a burst located at a proper 


location with respect to the earth's magnetic field could be 


detected from "horizon to horizon". 
This fact is illustrated 


in Figure 24. 
As one can see from this figure, the size of 


the area in which the EMP could be detected, for a given 


orientation with respect to the earth's magnetic field, 


depends on the height of the weapon burst. 


Calculations similar to that mentioned at the beginning 


of this section were performed for various gamma yields and 


heights of burst for the weapon. 
These calculations also 


ignored conductivity and considered only the gamma induced 


(ignoring x-ray induced) 
EMF. 
The results of these are 


presented on the following pages in graph form. 


69. \ 


GNE/PH/74-7 


E 
o 


\ 


C 
OBSERVERS 
B 


/ 


WEAPON BURST' 
I 


/~ 


// 
\ 
/ 


/ 
.; 


/ 


A 


Figure 23. 
Effect of the Earth's Magnetic Field' Orientation 
on the High Altitude EMP. 
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Figure 24. 
The Horizon-to-Horizon Nature of the EMP. 
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V. 
Effects of Air Conductivity 


In the previous chapters it was assumed that the con- 


ductivity of the air surrounding the radiating electron could 


be neglected. 
The purpose of this chapter is to study the 


validity of this assumption. 
This will be done by briefly 


studying the physics involved in air conductivity with re- 


spect to the 
m~chanisms relating to the creation of the 


conductivity and its subsequent buildup in time. 
There will 


also be a short section included to indicate the equations 


which must be solved in order to obtain the fields radiated 


within a conducting media. 


Formation and Buildup of Conductivity 


The basic cause of conductivity in the air is the 


ionization of the molecules of the air caused by the high 


energy Compton electrons. 
The result of this ionization is 


to release large numbers of low energy electrons 
(secondary 


electrons) into the media. 
Although these electrons are not 


energetic enough to add appreciably to the radiated fields, 


they do cause the environment to act as a conductor since 


they are free electrons. 
One can express the conductivity 


of the air in terms of the number density of these secondary 


electrons as follows 
(Ref 13, 
p 101): 


(V-I) 
a(R,t) = 
N 
(R,t)e 2 
sec 
m v 
(z) 
o c 


where a is the conductivity measured at a distance R from the 


weapon detonation at time t, 
N 
is the secondary electron 
sec 
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density, 
e is the magnitude of the charge on a single electron, 


mo is the mass of an electron, and 
Vc is the frequency of 


collisions at an altitude z above the ground. 
It should be 


mentioned that, in general, the conductivity depends on the 


density of negative and positive ions as well as free elec- 


trons. 
However, in this case the contribution of the ions is 


negligible as compared to that of the electrons. 
The reason 


for this is that the conductivity is proportional to the sum 


of the inverses of the mass of the electron, positive 
io~s, 


and negative ions. 
Since the ions are orders of magnitude 


more massive than the electron, the terms involving the 


inverse of their mass will be negligible compared to the 


term involving the inverse of the electron mass. 


If one assumes that v 
is constant with respect to time, 
c 
then the only time dependence in 
(V-I) is that in the Nsec 
term. 
There are several mechanisms which affect the time 


dependency of the density of the secondary electrons. 
One 


such mechanism is, of course, the rate at which energy is 


deposited by the Compton electrons to ionize the air, there- 


by creating the secondary electrons. 
However, the chemistry 


of the air is such that there are also mechanisms whereby 


these secondary electrons may be removed. 
The mechanisms 


include the attachment of a free secondary electron to an 


02 atom to form an 02- ion and the recombination of an 


02+ ion and .a N2+ ion with a free electron. 
Since the main 


purpose of this chapter is to test the validity of ignoring 


the conductivity, the removal mechanisms for secondary elec- 


trons will be ignored in order to arrive at a "worst possible 
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case~ for conductivity. 
Thus, 
a lower limit for the time 


during which conductivity can be ignored will be developed. 


tn order to relate the number of secondary electrons 


created to a given amount of energy deposited by the Compton 


electrons, one can define an average ionization energy W 


which is defined as the ratio of the energy lost by the 


Compton electron to the total ionization produced by it. 
That 


is, W represents the average energy required to be deposited 


in order to create a secondary eiectron. 
An accepted value 


for W in air is 34 eV 
(Ref 19, 
p 5-6). 
If a Compton electron 


was created at some time T and if 
~DEP(tC) represents the 


amount of energy which this electron has deposited at some 


time t 
after its creation, then the number of secondary 
c 
electrons at time t is given by 


(V-2) 


~DEP(t-T) 


Nsec(t) = ---W~---- 


It is obvious from this argument that the number density of 


secondary electrons at a distance R from the burst at a given 


time t 
will be dependent upon the number density of the Compton 


electrons at this point in space and time. 
That is, at any 


time t after detonation the number of secondary electrons at 


a certain R will depend on the number of all the Compton 


electrons which have been created previous to t 
and the amount 


of energy which each has deposited between the time of its 
dN 
,(R, T) 
t ' 
d t 
If 
prl 
h 
f 
h 
crea lon an. 
-dt 
represents t 
e rate 
0 
c ange 


of the primary electron number density at Rand T, then one 


can express the secondary electron number density at any time 
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t 
as follows: 


(v-3) 


t 
dN 
. (R, t) 
N 
(R , t) = f 
{--,p,,-r:;-:l 
_ 
sec 
dt 
o 


. SDEP (R, t-T) 
}{ 
W 
}dT 


where an R dependence has been included in sDEP to account 


for any variation with altitude. 


In order to obtain numerical values for 
(V-I) 
one needs 


an expression for v 
(z). 
The expression for v 
can be shown 
c 
c 
to vary with altitude primarily as the density of the air 


varies with altitude. 
Thus, if one assumes that the atmos- 
• 
pheric density varies exponentially with altitude then one 


can write: 


(V-4 ) 
-z/z 
v 
(z) = v e 
0 
c 
0 
where Vo is the collision frequency at sea level and Zo is 


the scale factor introduced in Chapter II. 
Using geometry 


indicated in Chapter II one can write the following for 


(V-4) 


v 
(z) = v e-(HOB-R)/Zo 
c 
. 
0 


Using 
(V-I), 
(V-3), and 
(V-4) calculations were made for 


the conductivity as a function of time at various altitudes 


and for weapons with different amounts of gamma yields. 
A 


functional dependence for sDEP(R,t) was obtained at each 


altitude from data such as that indicated in Table 3. 
The 


data was "fit" using a least squares technique and the re- 


suIting functional form was used in the integral in 
(V-3). 


The results of these calculations for an altitude of 40 Km 


and varying gamma yields are presented in Table 4 on the 


following page. 
To illustrate the variation in altitude, 
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the conductivity at a-given instant of time is plotted 


-as a function of altitude in Figure 
28~ 


Relative SignIficance of Conductivity 


To gain a deeper understanding for the significance of 


the values listed in Table 5, it is necessary to consider 


Maxwell's equations as applied to a conducting media. 
With 


regard to the equation which involves the current density, 


(V-5) 
- 
1 
V x B = 
~ 
J 
+ -- 
V 
c 2 


(V-6) 


one can express the current density as the sum of the Compton 


current density and the.current density of the secondary 


electrons as follows: 


J = J 
+ J 
comp 
sec 


In the previous chapters, it has been assumed that 
J 
~O. 
sec 
If one assumes that Ohm's law applies to the secondary elec- 


tron current density, 
namely that 


Jsec = crE 


then 
(V-5) 
can be written as: 


- 
1 
aE 
V x 
B - 
II 
J 
+ llv0E + --2 
'"I t 
v comp 
c 
a 


Using 
(V-6) 
and the remainder of Maxwell's equations as 


___.indicated in Chapter III and Appendix J3, one can obtain 


wave equations for the vector potential A and the scalar 


potential V such that: 18 


(V-7) 
{ 
2 
a 
V -ll 
0 -- - 
v 
at 


'. 
,-].1 
J 
I 
1 
a2 
A 
A 
A 
v comp 
-- -} I I=LI 
1= 
c 2 at 2 
V. 
v 
_1:- 
p 


Ev 
C 


18 
. 
h' 
. 
h 
d 
Assum1ng t at 
0 
1S constant W1t 
respect to space an 
time. 
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Of ALTITUDE 
fOR A fiXED 
TIME. 
· 
o 


0 
N 
. .. 
· 
0 


Q 
I 
Ow 


~...-1 
)IE · 
0 


,...... 
(J) 
:::s:::: 
E'N 
...-1 
....., . 
>-0 
t-- 
t-4 
:> 
.. 


t-4 
t-- CO 
uC: 
=:J 0 
a:z: 
0u 


~ 
0· 
0 


oo 


64.00 
56.00 
32 • 0 0 
40 • 0 0 
4-B• 00 
ALTITUOE(KM) 
24.00 


• ..Jr-------r--c::::::...----.,-----,.-----,.---~::=>-..,..._---___, 
'?6.00 


80 
. \ 


GNE/PH/74-7 


, 
If one applies a Fourier transform with respect to time to 


the left hand side of 
(V-7) 
the time derivatives can be done 


explicitly to yield the following for the operator L enclosed 


in the brackets: 


(V-8) 


fields. 
If one divides 
(V-8) 
by 


where LFT indicates that the Fourier transform has been applied, 


i=!=I; 
and 
W is the transformation variable. 
It can be shown 


that w corresponds physically to the frequency of the radiated 
(iw)2 
---z- the results are: 
c 


where 
c2=1/~ 
E 
has been used. 
Thus, one can obtain an esti- 
v v 
mation of the relative importance of the term involving 


conductivity by comparing the magnitude of the term ~ 
EVW 


to 1. 
Since the actual 
field~ which are radiated form a pulse 


which consists of many frequencies, it is necessary to arrive 


at some "characteristic frequency" of the pulse in order to 


obtain numbers for the ~ term for various values of o. 
A 


EvW 


crude approximation which is consistent with the Uncertainty 


Principle is to let 
W be equal to the inverse time width of 


the radiated pulse. 
In Chapter IV it was shown that the 


-8 
7 
pulse was about 5xlO 
seconds long. 
Thus, 
w=2xlO 
inverse 


seconds. 
The results of such calculations are indicated in 


the last column of Table 4. 
If the number .01 is arbitrarily 


set as the level at which ~ is to be considered significant, 


EvW 


then one can see from Table 4 that the conductivity can be 
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ignored for a gamma yield of one kiloton during the times 


which the high frequency EMP is being radiated. 
However, as 


the gamma yield is increased to five and then ten kilotons, 


the conductivity does become significant at relatively early 


times. 
Therefore, it is fairly safe to ignore the conduc- 


tivity for weapons with relatively low gamma yields (2 1 KT). 


One should be aware that this last statement is based on 


ignoring .01 as compared to 1. 


Field Equations for a Conducting Media 
• 


Equation 
(V-7) 
can be solved formalistically as was 


indicated in Chapter III using the Green'sfunction approach. 


That is, 


(V-9) I 
A(r,t) I 
V (r, t) 


where G(r,tlrl,t') is the Green's function and which must 


satisfy the following: 


(V-IO) 
{ 
2 
d 
I 
d 2 
} 
- 
1- 
- - 
V -~ 
cr -- - -- --- 
G(r,t rl,t'} 
= o(r-rl)o(t-t ' } 
a 
at 
c 2 at 2 


Equation 
(V-10) 
can be solved using Laplace transform 


techniques to yield 
(see Appendix C) 


(V-II) 


where Illxl is the modified Bessel function with argument x and 


ulxl is the unit step function with argument x 
(Ulxl=o for 
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x < 0 and ulxl = I for x 
> 0). 
The reader should note that 


the above expression 'was derived assuming that a was con- 


stant with respect to time which, 
from Table 4, is seen not 


to be the case. 
Yet, even with this simplifying assumption, 


the author was unable to successfully integrate 
(V-9) using 


(V-II). 
This result was included to indicate to the reader 


the difficulties involved in applying the microscopic model 


of the high frequency EMP to conducting media . 


.. 


83 \ 


GNE/PH/74-7 
i 


VI. 
Final Remarks 


T~e purpose of this thesis has been to present a model 


for the high altitude, high frequency EMP which could be used 


as a pedagogical tool. 
As such, the emphasis has been placed 


on developing and outlining the theory rather than obtaining 


state-of-the-art numerical values for threat levels of the 


EMP. 
However, this author feels that this model has great 


potential for application in this area. 
The reason for this 


is that this model allows one to see the effect of changing 


a greater 
n~ber of physical parameters than does the macro- 


scopic approach which was briefly described in Chapter I. 


For this reason, the following recommendations have been 


made to anyone interested in following up on the work pre- 


sented in this thesis: 


1) 
An attempt should be made to solve the field 
equations for a conducting media as indicated in 
Chapter V. 
The approach to this might involve: 
a) 
finding a numerical solution to the 
integral involving Green's function given 
in Chapter V in equations 
(V-9) 
and 
(V-II) 
or 
b) 
approaching the problem using a technique 
other than that presented in Chapter V. 


2) 
If one is able to solve Recommendation #1, then it 
will become necessary to make more precise calcu- 
lations for conductivity than were included in 
Chapter V. 
One should include the various loss 
mechanisms for the secondary electrons such as 
attachment and recombination. 


3) 
More study should be made regarding the cross- 
sectional line-of-sight area A(R). 
This could 
include sensitivity studies with respect to 
increasing or decreasing A(R) 
and the possibility 
of including the n dependence indicated in (III-40). 
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Appendix A 


Gamma Ray Interaction with Matter 


Attenuation Law 


The probability of .a photon interacting with the atoms 


of its environment can be expressed in terms of either the 


microscopic cross section or the macroscopic cross section of 


the material which constitutes the environment. 
The micro- 


scopic cross section 
(ac ) is defined as the effective cross 


sectional area which an atom presents to a gamma photon· for 


interaction. 
The macroscopic cross section 
(~) is a more 


useful quantity than the microscopic cross section than 


driving a law for the attenuation of gamma radiation. 
The 


macroscopic cross section is defined as the probability of 


reaction per unit length of photon path and is related to the 


microscopic cross section by the relation 


{AI} 
~ = Nac 
{Ref 1, 
P l-6l} 


where N is the atom density of the environment. 
Since a 
has 
c 


{A2} 


units of area and N has units of inverse 
vol~~e, then 
~ must 


have units of inverse length. 
It is obvious from 
(AI) that 


~ will depend upon the physical state of the environment at 


any given time, whereas a 
is a fundamental quantity which is 
c 
constant for a given material regardless of the physical state. 


The quantity N can be found from the following relation 


~aP 


N = AW 


where Na is Avogadro's number, 
p is the material density, and 


AW is the atomic weight of the matter making up the environment. 


For the purpose of this study an attenuation law with 
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respect to a spherical geometry 
(see Figure 29) will be 


developed. 
As shown in Figure 29, the source shall be con- 


sidered a point source such that photons emanate from it 


equally in all directions. 
Thus the only variation in the 


photon intensity (photons/unit area) is with respect to the 


distance from the source r. 
If 
~ 
represents the number of 
o 
photons which originate from the source and if one assumes 


that there is no removal mechanism for the photons in the 


transport media then the number of photons crossing a unit 


area at a distance r 
from the source (intensity - 
~(r» 
is 


given by 


(A3) 
~(r) = 


since 4rrr 2 represents the surface area of a sphere of radius 


r. 
If a material of cross section 
~ is placed around the 


source as a shell with an inner radius of r 
and thickness dr 


(see Figure 29A) 
then the intensity of unreacted photons 


(virgin intensity) at r 
+ dr will 
b~ 


(A4) 
~(r+dr) = 
~(r) 
+ 
d~(r) 
(see Figure 29B) 


In the previous equation 
d~ represents the change in the 


virgin photon intensity across dr because of interaction 


with the material shell. 
Since 
~ represents the probability 


of an interaction per unit photon path length, this change 


d~ must satisfy the relation 


(A5) 
d~(r) = 
-~(r)~dr 


That is, the number of photons which interact in the thick- 


ness dr is equai to the total probability of photon interaction 


in the shell 
(~dr) multiplied by the number of photons "entering" 
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Fig~:r'e 29. 
Spherical Geometry for Attenuation Law. 


the shell. 
The equation 
(AS) 
can be rewritten such that 


(A6) 
dq,(r) 
= -lldr 
q,(r) 


If one pictures the gamma source as completely surrounded by an 


. 
.... 
infinite numbe'r of shells of infinitesimal thicknesses such that 


the inner most shell has an 
inner·~adius of zero and the outer- 


most shell has an outer radius of R then the virgin flux at R . 


can be obtained by integrating the left hand side of 
(A6) 


q,o 
between 
4nR L and q,(R) 
and the right hand side between 0 and 


R. 
That is 


q, (I~) 
R 


f 
dq,/q, = 
f 
lldr{assume ll=ll(r)} 


q, 
/4nR~ 
0 
o 
" 
or 


{q,(R~4nR2 
'R 
In 
} = 
f 
lld;r 


o 
' 
0 
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or equivalently 
R 


"':'1 lldr 


o 
= e 


or solving for 
q> (R) 
R 


q> 
-J lldr 


(A7) 
q>(R} 
0 
0 
= --e 
41TR 2 


The last expression above represents the attenuation law 


which is used in this thesis for the virgin intensity of 


gamma photons. 
• 


compton Effect 
(Ref 16, 
p 432) 


The primary gamma photon interaction with respect to the 


high altitude EMP is the Compton effect. 
The analysis of this 


effect is performed by applying the energy and momentum con- 


servation laws to a collision between an electromagnetic 


quantum 
(gamma photon) 
and an electron. 


of the incident gamma, 
E ' 
Y 


m c 2 the rest energy of the 
o 


If E 
is the energy 
y 


the energy of the scattered gamma, 


electron, and W the total energy 
e 
of the recoil electron 
(see Figure 6, Chapter II) 
then the 


energy conservation law can be expressed as 


Ev 
+ m c 2 = E ' + W 
lOy 
e 


Using the substitutions indicated in Figure 
6 in Chapter II 


one obtains for the above 


(A8) 
hv + m c 2 = hv' 
+ W 
o 
e 


Similarly, the momentum conservation laws for this situation 


can be written as 


(A9) 
hv 
c 
= hv' 
cos e + p 
C 
e cos 
~ 
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and 


(AID) 
hv' 
D = 
sin e - P 
sin 
~ 
c 
e 


where the variables used above are those indicated in 


Figure 6. 
The variable 
~ can be eliminated by rearranging 


(A9) 
and 
(AID), squaring both sides and adding. 
The results 


of this are 


= P 
2 
e 


If one multiplies the above by c 2 and then rearranges terms, 


the equation becomes 


(All) 


Equation 
(A8) 
can be solved for W 2 whereupon one-finds 
e 


(A12) 


If one subtracts equation 
(All) 
from 
(A12) 
and eliminates We 
and P 
through the relativistic relation 
e 
W 2 
_ P 
2 C 2 = m 
2 C 4 
e 
e 
0 
then one obtains 


(A13) 
2h(v-v')m c 2 = 2hvv' (I-cos 8) 
o 
Equation 
(A13) 
can then be easily solved for hv' with the 


result that 


hv' = 


or 


E ' = 
y 


{I + 


{I + 


hv 
hv~ 
o 


EY 
EY 
iiiC' 
o 


(1 - 
cos 8)} 


(1 - cos 8)} 
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Appendix B 


Derivation of E 'and B for an Accelerated Electron 


The Wave Equation 
(Ref 4) 


Maxwell's equations are the starting point for finding 


relations for E and E for a given source. 
These equations 


are 


(Bl) 


(B2) 


(B3) 


(B4) 


v x E(r,t) 
= 
dE(r,t) 
dt 


V x B(r,t) 
= 
II J{r,t) 
+ llvE:v 


dE{r,t) 
v 
~t 


VoE(r,t) 
1 
p (r,t) 
= 


E:v 
c 


VoB(r,t) = 0 


where the variables are as indicated in Chapter III. 
Hence- 


forth it shall be assumed that the conductivity of the media 


surrounding the' current source is zero. 
Because 
(B4) is 


always true one can say that 


(BS) 
B = V x A 


where A is an arbitrary vector at this point. 
Substituting 


CB5) 
into 
(BI) 
and assuming that A is such that the curl 


operator and time derivative operator can be interchanged, 


one obtains the following: 


(B6) 
V x {E(r,t) + dA{r,t)} = 0 
at 


Since for any scalar V it is true that VxVV=O, one can make 


the following assertion without loss of generality: 


- 
dA 
E + - 
= - 
VV 
at 


or 
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(B7) 
aA 
E = - at - 
VV 


where V is an arbitrary scalar. 
If one substitutes 
(B7) 


into 
(B3) 
and again assumes that A is such that the spatial 


derivative 
(Vo) 
can be 
in~erchanged with the time derivative 


then the following relation is found: 


(B8 ) 
2 
- 
a 
- - 
v V(r,t)+at {V·A(r,t)} = - 
p (r, t) 
c 


Also, if one substitutes 
(B7) 
and 
(BS) 
into 
(B2) 
and uses 


the identities Vx(VxA)=V(V·A)-V 2A and 
~ 
E =~ then one 
v 
v 
c 
obtains 


(B9 ) 


At this point there is still a certain "arbitrariness" with 


respect to A. 
Since any vector is uniquely defined when both 


its curl and divergence are specified, the arbitrariness in 


A can be removed by specifying. its divergence or the gauge. 


By studying 
(B8) 
and 
(B9) 
one observes that a convenient 


gauge for this problem is the Lorentz gauge or 


- 
I 
av 
(BID) 
VoA = - ex 
3t 


substituting 
(BID) 
into 
(B8) 
and 
(B9) 
one finds the following 


relations for the wave equations of A and V: 


(BII) 


and 


p (r,t) 
c 
' 


(BI2) 
I 
-2 
- 
~ J(r,t) 
v 


The Green's 
Functio~ (Ref 4) 


Solutions to 
(BII) 
and 
(B12) 
can be written formal- 


istically using the Green's function approa?h such that the 
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19 
following is true: 


(B13) . 
A(r,t) 


V (r, t) 


-llVJ(r' ,t') 
_L- 
p (i' 
t') 


E: 
C 
' 
V 


(B14) 


where G(r,tlr' ,tIl is known as the Green's function and 


which must satisfy the following equation 


1 
a2G(r,t r' ,tIl 
V2G(r,tlr',t) 
- cY 
at~ 
= o(r-r') 
oCt-tIl 


One should note that the derivative operators 
(V 2 and 
~~2) are 


both with 
re~pect to the unprimed coordinates. 
To solve 
1B14) 


one can use the Laplace transform technique to remove the time 


derivative from the problem. 
Thus, if one lets A denote the 


Laplace transform operator with respect to time then the Laplace 


transform of G becomes 


(B15) 


_ 
00 


G(r,slr' ,t')=A{G}=fdte-st G(r,tli;t') 
o 
. 
th 
f 
L 
1 
t 
f 
f 
h 
th d' 
. 
USlng 
e property 
or 
ap ace 
rans orms 
0 
t 
e p 
erlvatlve 


of a function F(t) which is as follows 
(Ref 8, 
p 490) 


n-l 
A{F (p) (t)} = sPf (s) 
- I 
sp-l-n F (n) 
(+0) 
n=O 


where f(s)=A{F(t)}, one can then write for the Laplace trans- 


form of the second time derivative of G the following: 


A{a2G(ra~!r',tIl 
} = 
s2~(r,slr' ,tIl 
- sG(r,olr' ,tIl 


_ aG (r, t Ir' ,t ') I 


at 
t=O 


If one demands that the Green's function and its first time 


19This is just a short hand notation to avoid writing 
both equations. 
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derivative satisfy homogeneous boundary conditions at t=O then 


the previous equation becomes: 


A{Cl 2G(r,t/r' ,t') 
} = s2G(r,slr' ,t') 


Clt 2 


One can find the Laplace transform of the right hand side of 


equation 
(BI4) directly by solving the equation: 


00 
A{o(r-r')o(t-t')} = f dte-st oCt-tIl 
o 


which, if one applies the properties of the Dirac delta 


function, 
becomes 


(BI7) 
A{o(r-r')o(t-t')} = o(r-r') e-st ' 
Using equations 
(BI5), 
(BI6) 
and 
(BI7) 
one can see that if 


the Laplace operator is applied to both sides of equation 


(BI4) that this equation becomes 


(BIB) 


where the r subscript to 


S2 
- 
- 
- 
-st' 
- -- G(r,slr',t') 
= oCr-rIle 
c 2 


the V2 operator is used to emphasize 
r 
that the derivative is with respect to the unprimed r. 
By 


translating the origin from 0 to r' 
and letting R=r-r' 
one 


can write 
(B18) 
as 


(BI9) 


where the argument of G has been dropped for brevity's sake and 


where it has been assumed that v: 
+~ 
v~, (this is easily shown). 


If one multiplies 
(BI9) 
by est' 
on both sides and lets G=G est' 


then this 


(B20) 


equati<;m 


~ 
S2:::: 


G - -- G 
c 2 


can be written as follows: 


= 0 (R) 


Since oCR) represents essentiaLly a point source at R=O, one 


can argue that there is spherical syrrmetry in the problem and 
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the only spatial dependence which G has is with respect to the 


distance from the source and not on direction. 
That is, 


G= 3(IRI) = G(R) 
- 
Thus for 
V2 G(R) 
one can write the following 
R 
-- 
V2 
G 
R 


- 


= .!- ~ 
(R 2 dG) 
R2 dR 
dR 


Using the above and letting G=g/R one can rewrite 
(B20) as 


1 d 2 g 
_ 
S2 
g- = 0 (R) 
R dR2 
c 2 
R 


or 


(B2l) 


If R>O then o(R)=O and the right hand side of 
(B20) is zero; 


if R=O then the integral of the right hand side is zero. 
Thus, 


without loss of generality, one can solve the following instead 


of 
(B21): 


(B22) 


Equation 
(B22) 
has the following general solution: 


Since it is desired that the fields go to zero as R approaches 


infinity, then C2 must be zero in the above. 
Also, in order 


that the final solutions would reduce to the electrostatic 
-1 
20 
solution for V(r,t) it is seen that Cl must equal 4TI' 
Thus 


20This is stated without proof. 
The interested reader 
can carry Cl 
through the derivation to find V and then match 
the solution to the electrostatic one to find Cl . 
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one has for g(R) that 


1 
-sR/c 


= - 
41T e 


- 
then for G, 


-G = 


and for G, 


(B23) 


-sR/c 
g(R) = -e 
R 
41TR 


-st' 
:::; 
= e 
G = 


A general 
id~ntity concerning Laplace transforms is that· 


(B24) 


where A-I denotes the inverse transform. 
Therefore, if one 


lets T=t,+Jr-r' '- and applies equation 
(B24) to 
(B23) 
the 
c 
result is 


(B25) 


-s (t,+Jr-r' '-) 


- 
1 
A-I{e 
C} 
41Tlr-r'l 


= - 
1 
6{t-(t' + Jr-r' L)} 
41T Ir-r' I 
c 


Since the Dirac delta function is an even function and since 


written as 


(B26) 


Going back to equation 
(B13) 
one then has for A(r,t) 
and 
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V{r,t) 
the following relation: 


(B27) . 
A (r, t) 


V (r, t) 


llVJ (r' , t ' ) 


1 
pc(r' ,t') 


£v 


E and B Solution 


Before one can find a solution for E(r,t) 
and V(r,t) 


from 
(B27) 
one must have an analytical expression for J(r' ,t') 


and pc(r',t'). 
For an electron with charge e and velocity 


vet) at the point x(t) 
the current and charge density are· 


given by the following: 


and 


where S(t')=v(t')/c. 
Thus one can write for the vector and 


scalar potential the following: 


V (r, t) 


]JvcB(t') 


(){r'-x(t')} 
1 


£v 


using the properties of the Dirac delta function the integration 


over the spatial coordinates (r') 
can be immediately performed 


to yield: 


(B28) 


A(r,t) 
o{t,+R(t') -t} ]Jv.cB(t') 
=~fdt' 
c 
V(r,t) 
4'lT 
R(t') 
1 
£v 


where R(t') 
= Ir-x(t') I 


One should observe that the only spatial dependence in 


(B28) is in the term R. 
Further, since the forms of E and B 


given by equation 
(BS) 
and 
(B7) require the use of the spatial 
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operator V 
it is convenient to express this operator in terms 
r 
of its effect of R rather than r. 
lising the Einstein 
su~nation 


, 
21 
notation 
the vector R can be written as 


R = x.e. 
1 
1 
(i = 1, 2, 
3) 


where the x. represent the coordinates of R in the cartesian 
1 
basis 


(B29) 


e .. 
1 


Vr 


One can also represent the Vr operator by 


" 
a 
-+ e. 
- 
1 
ax.1 
so that for some function f(R) 


v f' (R) 
= 
~. af (R) 
r 
1 
ax.1 


or, using the chain rule 


" 
af (R) 
aR 
V f (R) = e 
r 
i 
aR 
ax.1 


" 
aR 
af(R) 
Vrf(R), = e. 
1 
ax. 
aR 
1 
or, using 
(B29 ) 


= f{R(x.)} one can write that 
1 


• 
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pointing from the source to the observer. 
Thus one can write, 


using 
(B30) 
and 
(B3l), that 


(B32) 


Using the relation for E given in 
(B7) 
as applied to 
(B28) 


and letting 
~ 
be as given in 
(B32) 
one finds the following 
r 


= :'IT fdt' 
Co: 
R~ (t') 
o{t' 
+ 
v 


for E (r, t) 


E(r,t) 
R(t') 
_ t} 
c 


A 
a(t' 
+ R(t') 
-t) 
n 
{o' (t,+R(t') 
-t.) 
-=c:---__• 
} 
s R(t') 
c 
aR 
v 


B(t')~vc 
R(t') 
a(t'+R(t')/c-t 
{o' (t'+ 
-t) 
}) 
R(t') 
c 
at 


where 0' (x) denotes the first derivative of the Dirac delta 


function with respect to its argument. 
After terms are re- 


arranged and the indicated derivatives performed the above 


expression for E(r,t) 
can be written as 


(1333) 
··-E(- t) 
= _s_ fdt' (n 
~{t' 
+ R(t') -t} 
. r, 
4 
v 
'IT Sv 
R2 (t ' ) 
c 


+ 
1 
(S(t')-~) o'{t' 
+ R(t') -t}) 
R(t')c 
c 


where the fact that 
~ 
E 
= 1 
has been used. 
In order to 
vv 
2" 


evaluate the integral in 
(B33) it will be divided into the 


following two parts: 


(B34) 
I 
= fdt' 
1 
o{t' 
+ R(t') 
-t}) 
c 


(B35) 
I 2 = fdt' 
R(~')C 
(S(t')-~) {o' (t' 
+ R~t') -t)} 


It should be noted that the argument of the delta function 


in II (andI 2 ) contains a function of the_ variable of integration 
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rather than simply the variable itself. 
In order to evaluate 


this particular situation the following general integral will 


be considered: 


I = fg(t l ) o{f(t l ) -c}dt' 


where g and f are arbitrary functions of x each of which has 


an inverse22 and first derivative and where c is some constant 


with respect to the integration variable. 
If one lets f(t')=u 


.. 


du = df 
dt' 
4 dt' 
= 
dt' 


and I 
becomes 


du 
df/dt' 


(f- 1 ) 
I 
= f 
~fldt' 
o(u-c)du 


or using the delta function properties 


or 


(B36) 
I = fg(x)o{f(tIY-c}dt l 
= r g(ti 
l ) 
i 
df/dt' f(t')=c 


where the t. I are those values of t l which satisfy f(tl)=c. 
1 
Applying 
(B36) to 
(B34) 
one can see first that f(tl)=tl+lR(t') 
c 
which implies the following: 


(B37) 


If one represents R by 
(R'R)1/2 then ~ ~~l can be written as: 


1 dR 
c dt' 


A_R 
or using the definition of R and letting n=R 


22For example if f(x)=u, 
then x=f·· 1 (u) 
where f- 1 is the 
inverse function of f. 
loi \ 
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(B38 ) 
1 dR 
c dt' 


A 
1 d 
- 
- 
A 
V (t' ) 
= n o- 
- 
{r-x(t') }=-n·o--,---,-= 
c dt' 
c 
noS(t')" 


Thus one has for .(B37) 


(B39) 


Using 
(B38) 
and the relation defined in 
(B36) 
one finds that 


(B34) 
can be written as follows 


or 


I 
= 
I 


A 
n 
t' + R(t') 
= t 
c 


(B40) 


A 


It'=t _ R(t') 
= I ~ I 
C 
--2 
ret 
KR 


The brackets with the subscript "ret" means that the quantity 


which they enclose is to be evaluated at the retarded time, 


t'=t-Ro 
c 
The expression 
(B35) 
and I 2 can be evaluated similarly to 


II if one makes the following substitutions: 


f (t') 
= t' + R(t') 
_ u 
+ t' 
= f- 1 (u) 
C 


(B41) 
dt' 
du 
du 
= df/dt' 
- 
K 
( as before) 


Then one has for I 2 that 


I 
= I 
J g (u) 
~'( 
t) d 
2 
c 
df/dt' 
v 
u- 
u 


The above can be integrated by parts on the derivative of the 


delta function to yield: 


I 
= _ ·1 Jo(u-t) 
d 
{g(u)} dt' 
du 
2 
c 
dt' 
K 
du . 


This can then be integrated in the same manner as II to yield 
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I, 


1 (dt' 
d 
I 2 = - 
c 
du 
dt' 
{g(U)}\ 
K 
u=t 


or using the substitutions indicated in 
(B40) 


(B42) 
I 
= - IL 
2 
KC 


A 


~t ' 
(_~;n) Iret 


The equation 
(B33) 
for E(r,t) 
can be written in a closed 


form if one employs 
(B40) 
and 
(B4l). 
The results of this are 


E (r, t) 
e 
= 4'ITEv 
+ 1 
d 
KC dt' 


1 dn 
To perform the differentiation - 
c dt' 
one can use the following 


identity presented here without proof or justification23 


A' 


(B43) 
1 dn 
= 
c dt' 


Using the substitution indicated by 
(B43) 
and performing the 


necessary algebra, one finds the following relation for 


E(r,t): 


(B44) 


A 
A 


E(r,t)=_e_/ 
n 
+n 
L.(L)_~ 
1 
d 
4'ITE 
2R2 
CK dt' 
KR 
2 
2 
-CK dt' 
v 
K 
K R 
ret 


If one had performed similar calculations to those performed 


up to this point on B(r,t) where B=VxA then the following 


relation would have been found for B(r,t): 


ec~ 
- 
- 
(B45) 
B(r,t) 
= ----4v {(_s_ + 1 
d 
{S}) 
x 
~} 
'IT 
CK dt' 
KR 
K 2 R2 


One can easily show using 
(B44) 
and vector algebra that the 


following relation exist between E and 
B~ 


(B46) 
_ 
_ 
1 
A 


B(r,t) 
= - 
{n x E(r,t)} 
. 
C 


23 
'f 
To verl y 
r=x(t')+R(t') 
and 
the definition of 


this relation one should realize that 
dr 
that for a fixed observer dt-O. 
Also use 


~=R(R.R)-1!2. 
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In order to obtain a more explicit form forE(r,t) 


than that given by 
(B44) 
the following relation is noted: 


A 


(B47) 
1 d 
R 
A 
~ 
R- dn 
1 dR 
A 
- 
c dt' (KR)= -c(noS)-cSOdt' 
+ c dt' (I-noS) 


where 


(B48) 


Equation 
(B47) 
can be rewritten using the substitutions 


indicated by 
(B38) 
and 
(B43) 
such that 


(B49) 
1 dAR A ~ 
- --- 
(KR) 
= S2 - 
(Son) 
- 
c(noS) 
c dt' 


Therefore, if one uses the substitutions given by. equations 


(B48) 
and 
(B49) 
and rearranges the terms using vector algebra 


then equation 
(B44) 
can be written as follows: 
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Appendix-C 


Derivation of the Green's Function for Conducting Media 


If one assumes that a conducting media obeys Ohm's law 


and that the conductivity is uniform throughout space and 


constant with respect to time then the following wave equation 


for the vector and scalar potential can be derived 
(see 


Chaptet V) 


(el) 
{ 
2 
a 
If -ll 
(J -- 
v 
at 


A(r,t) 


V (r,t) 
= 


-llVJ(r,t) 


1 
- 
-£ pc(r,t) 
v 


(C2) 


The Green's function associated with this must satisfy the 


following equation 


2 
a 
1 
a2 
{If -ll 
(J- -- -- 
} G(r,tlr' ,t') = o(r-r')o(t-t') 


v at 
c2 at 2 


One can solve 
(C2) 
in exactly the same manner as indicated in 


the second section of Appendix'B using Laplace transform tech- 


niques. 
Doing this one finds the following for the Laplace 


tra.nsform of G(r, t Ir' , t' ) : 


(C3) 


.... 
G(r,slr',t') 
= 
-st' 
e 


It will be necessary to rearrange some of the terms and use 


some transform identities in order to get 
(C3) 
into a form 


which can be found in a table of Laplace transforms. 
The 


following identity shall be employed 
(Ref 8, 
p 491): 
-st' 
(C4) 
A- 1 {f(s)e 
} = F(t-t') 


where 


A- 1 {f(s)} = F(t) 
and F(t) = 0 for t<O 
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Therefore it is only necessary to find the inverse transform 


of the following for 
(C3): 


r 
c 
.; s (~ + s) 
(CS) 
g(s) = e 
t.v 


If one lets K=~ then one can write that 
t.v 


s(s + ~) 
= s(s + K) 
t.v 


or 


s(s + K) 


or 


(C6) 
s(s + K) 
= 
(s + 
~) 2_(~) 2 


Using 
(C6) 
one can then write the following for 
(CS): 


~ ;I(s + ~)2 _ 
(~)2 


(C7) 
g(s) 
= e 


If F(t} = A-1{h(s}} then the following identity is true 


(Ref a, 
p 491): 


(ca) 
A-1{h(s - 
b}} = ebt f(t) 


where b is a constant with respect to the transform. 
Thus, 


if one lets h(s} 
be the following: 


- 
~ 
/ 
S2 
_ 
(~) 2 
h(s} = e 
then from 
(C7) 
r 
c 
h(s-b) 
= g(s) 
= e 


K 
where, in this case, b=-2. 
Thus, using 
(Ca) 
one finds 


(C9) 


_!S.t 


A-1{g(s)} = A-1{h(s-b)} = e 
2 


_r.; S2 
_ 
(K)2 
1 
C 
2 
A- {e 
} 


The last term in 
(C9) 
can be written as follows by adding 
o to both sides: 
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r 
c 


(CIO 
-1 { 
A 
e 
rs 
-- 
c 


-e 
} 


-I{ 
+ A 
e 
} 


The inverse transforms of each term on the right hand side of 


(CIO) 
are common in tables of Laplace transforms. 
Using such 


a table the following is found: 


(Cll) 
-1 { 
A 
e 
} = 
ult-r , 
c 


r 
+ 0 (t--) 
c 


where Illxl is the modified Bessel's function with argument 


x and ulxl is the unit step function with argument x. 
Sub- 


= e 
-1 
A 
{g(s)} 


stituting 
(Cll) 
into 
(C9) 
one finds that: 


-~t 
2£v 


(C12) 


Finally if one applies 
(C12) with the identity given in 
(C4) 


to 
(C3) 
the result is as follows: 


x ult-t'-£l 
c 


'107 \ 


(C13) 


GNE/PH/74-7 


Vita 


John R. Lillis was born on 4 November 1950 in Frankfort, 


Kentucky. 
He was graduated from Madison High School in 


Richmond, Kentucky in 1968. 
He attended Virginia Polytechnical 


I~stitute, Georgia Institute of Technology, and received the 


degree of Bachelor of Science with a major in physics and a 


regular commission in the United States Air Force from the 


University of Kentucky in May, 1972. 
After graduation he was 


assigned to the Air Force Institute of Technology where he 


received a Master of Science degree in Nuclear Engineering 


in March, 1974. 
He is a member of Phi Beta Kappa and Tau 


Beta Pi academic honoraries. 


Permanent mailing address: 
1421 Fairlane Drive 
Richmond, Kentucky 40475 


This thesis was typed by Mrs. Laura Reams. 


108. \ 


